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ABSTRACT 


This  thesis  develops  the  mathematical  expressions 
to  evaluate  component  weight  and  space  for  all  signifi- 
cant components  within  the  containment  vessel  (secondary 
shielding)  of  a  Gas  Cooled  Reactor  (GCR)  closed  cycle 
gas  turbine  plant  for  ship  propulsion.   The  equations 
are  solved  for  a  range  of  operating  conditions  using 
helium  as  the  working  fluid.   Information  from  these 
equations  was  used  to  obtain  the  following: 

1.  Optimum  total  weight  of  components. 

2.  Curves  of  component  weight. 

J.  Curves  of  component  dimensions. 

The  resulting  data  support  the  following  conclu- 
sions i 

1.  Optimum  weight  occurs  at  a  pressure 
ratio  that  is  lower  than  the  value  that  gives  best 
efficiency. 

2.  A  heavy  weight  penalty  must  be  paid 
for  small  increases  in  efficiency, 

3.  The  regenerator  is  the  controlling  com- 
ponent in  a  cycle  of  this  type.   As  heat  transfer  re- 
quirements go  up  (regenerator  effectiveness  approaches 
100  %),  regenerator  weight  rises  sharply  and  limits 
appreciable  gains  in  efficiency. 

4-.  Rotating  machinery  length  will  have  a 
greater  influence  on  plant  design  than  the  weight  of 
this  component. 

The  curves  should  give  the  designer  a  rough 
idea  of  machinery  component  sizes  and  weights,  and 


li 


are  intended  for  use  in  the  development  and  design  of 
an  efficient,  nuclear  gas  turbine  power  plant. 
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NOMENCLATURE 

A         cross  section  or  heat  transfer  area  (ft2 ) 

C         absolute  velocity  (ft/sec) 
constant 

Cp        specific  heat  at  constant  pressure  (BTU/Lb  R) 

D         diameter  (ft) 

d         pitch  diameter  at  mid-blade 
diameter 

,         .  ,   -,  .   ,.    .     4-(cross  section  area)   /^4_\ 
de        hydraulic  diameter  =   Netted  perimeter   '  (ft) 

d         equivalent  outside  flow  diameter  per  tube  (ft) 

f         Fanning  friction  factor 
blade  taper  factor 

g         gravitational  constant 

H         surface  film  coefficient  of  heat  transfer  (BTU/ft2 -HR-°F) 
height  (ft) 

J         conversion  factor  (778  ft-Lb/BTU) 

K         ratio  of  specific  heats 
constant 

k  thermal  conductivity  ( BTU/ft -HR-°F) 

L  length  (ft) 

j£  blade  length 

M  mach  number 

m  mass  (Lb) 

N         number  of  tubes 

number  of  fuel  rods 
number  of  stages 

P         pressure  (psia) 

p         distance  between  tube  centers  (tube  pitch),  (ft) 
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Pr       Prandtl  number 

Q        heat  rate  (BTU/HR) 

q        constant,  function  of  tube  arrangement 

R        gas  constant  (ft-Lb/Lb°R) 
radius  (ft) 

r  pressure  ratio  (>l) 

Re  Reynolds  number 

S  entropy  (BTU/Lb°F) 

T  temperature  (°R  =  °F  +  460)  and  °F 

U        heat  transfer  coefficient  (BTU/ft2 -HR~°F) 
circumferential  velocity  on  pitch  diameter 


u 


exponent  =  j?— r 


V  velocity  (ft/sec) 
volume  (ft3 ) 

K— 1 
v        exponent  =   -g— 

W        work  (BTU/Lb) 
weight  (tons) 

w        flow  rate  (Lb/HR) 

Y  pressure  loss  factor 


vv 


£        stress  parameter  =   — —   (ft) 
angle  of  relative  velocity 

6        blade  stress  (p  si) 
r\  efficiency 

^x       regenerator  effectiveness 
-J\  angular  velocity  (rad/sec) 

ix 


t         ratio  of  compressor  isentropic  exit  temperature  to 
its  inlet  temperature 

^  density  (Lb/ft3)  or  (tons/ft3) 

V-  viscosity  (Lb/HR-ft) 

A.   T       log  mean  temperature  difference 


Subscripts 


A 

accumulator 

a 

axial 

area 

c 

compressor 
cold 

d 

diameter 

S 

gas 

HPC 

high  pressure  compressor 

HPT 

high  pressure  turbine 

h 

hot 

i 

inside 

ic 

intercooler 

LPC 

low  pressure  compressor 

LPT 

low  pressure  turbine 

JL 

length 

m 

mean 

N 

number  of  tubes 

0 

outside 

P 

piping 

pc 

precooler 

R 

reactor 

RG 

regenerator 

s 

surface 

stage 

t 

turbine 

V 

reactor  vessel 

w 

water 

weight 

cross  section 


State  Points 


1.  precooler  outlet 
LPC  inlet 

2.  IPC  exit 
intercooler  inlet 

2s o       isentropic  LPC  exit 

3-        intercooler  outlet 
HPC  inlet 

4.  HPC  exit 

cold  side  regenerator  inlet 

4s c       isentropic  HPC  exit 

5.  cold  side  regenerator  outlet 
reactor  inlet 

6.  reactor  outlet 
HPT  inlet 

7.  HPT  exit 
LPT  inlet 

7s.       isentropic  HPT  exit 

8.  LPT  exit 

hot  side  regenerator  inlet 

8s.       isentropic  LPT  exit 

9«        hot  side  regenerator  outlet 
precooler  inlet 
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I.  INTRODUCTION 

The  thermodynamic  potential  of  the  closed  cycle  gas 
turbine  power  plant  has  been  well  recognized  for  many 
years.   Early  attempts  to  effectively  utilize  this  poten- 
tial, however,  were  unsuccessful  because  of  poor  compon- 
ent efficiencies  and  the  lack  of  a  suitable  high  tempera- 
ture heat  source.   Technical  developments  of  high  tempera- 
ture gas  turbines  over  the  past  years  have  enabled  the 
designer  to  develop  components  with  the  necessary  thermal 
efficiency,  and  the  advent  of  the  nuclear  reactor  has 
created  the  type  of  heat  source  necessary  to  drive  the 
system. 

In  the  design  of  a  ship  propulsion  plant,  certain 
parameters  are  not  subject  to  the  designer's  control  or 
choice.   The  maximum  temperature  in  the  cycle  (T^)  may 
be  restricted  because  of  reactor  fuel  cladding  creep 
strength  limitations  or  the  restriction  imposed  by  the 
high  pressure  turbine  (HPT)  blading.   The  minimum  tempera- 
ture (T.)  of  the  working  fluid  is  likely  to  be  set  from 
the  sea  water  inlet  temperature  and  reasonable  size  and 
weight  coolers.   A  fixed  reactor  inlet  temperature  (Tj-) 
is  necessary  to  prevent  cyclic  thermal  stresses  from 
damaging  the  reactor  structural  material  and  moderator. 
In  a  gas  cooled  reactor,  the  inlet  temperature  is  put  at 
a  relatively  high  value  to  anneal  out  radiation  damage 
in  the  moderator.   Compressor  and  turbine  efficiencies 
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are  dependent  on  the  state  of  technical  development  of 
these  machines.   In  view  of  the  above,  T,-,  T-.,  Tj-,  and 

rotating  machinery  efficiencies  will  most  likely  be 
fixed  for  any  cycle  because  of  mechanical  design  or 
metallurgical  restraints  and  size  of  the  plant. 

The  designer  then  has  at  his  disposal  other  para- 
meters which  he  must  fix  so  as  to  give  the  most  economic, 
efficient  or  the  least  weight  and/or  space  propulsion 
plant o 

The  closed  cycle  nuclear  gas  turbine  plant  con- 
sidered in  this  work  consists  of  two  principal  sections: 
the  reactor,  which  is  located  within  the  primary  shield, 
and  the  propulsion  equipment  which  is  located  within 
the  secondary  shielding  or  containment  vessel  along 
with  the  reactor.   The  main  objective  of  this  thesis  was 
to  develop  a  weight  and  space  optimizing  procedure  for 
all  significant  components  located  within  the  contain- 
ment vessel.   The  machinery  components  considered  as 
significant  with  respect  to  the  variation  in  their  weights 
and  sizes  with  changes  in  the  cycle  parameters  ares  the 
regenerator,  coolers,  rotating  machinery,  reactor,  piping, 
and  accumulator.   The  developed  expressions  for  component 
weights  and  dimensions  in  terms  of  thermodynamic  proper- 
ties, gas  constants  and  flow  rate  enables  the  designer 
to  compare  the  relative  size  and  weight  of  each  component 
for  different  working  fluids. 
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The  second  thesis  objective  was  to  make  a  weight 
and  size  calculation  of  these  components  using  the 
proposed  helium-cooled,  graphite  moderated,  Maritime 
Gas-Cooled  Reactor  (MGCR)  power  plant  as  a  reference. 
Solution  of  the  developed  equations  for  a  range  of 
operating  conditions  provided  data  for  determining 
minimum  weight  and  space  for  the  regenerative,  inter- 
cooled  gas  turbine  cycle. 

No  gas  turbine  nuclear  power  plant  has  been  or 
is  presently  being  constructed  for  ship  propulsion. 
However,  with  the  technological  advances  being  made 
in  nuclear  reactors,  high  temperature  materials  and 
rotating  machinery,  it  is  to  be  expected  that  a  ship 
power  plant  of  this  type  will  someday  be  used.   The 
shortage  of  component  weight  and  space  data  for  such 
a  system  will  present  a  major  problem  to  the  ship  pre- 
liminary designer  who  must  develop  an  efficient  nuclear 
gas  turbine  plant.   The  developed  weight  and  space  equa- 
tions  along  with  the  system  study  of  a  particular  plant 
undertaken  in  this  thesis  should  assist  the  naval  archi- 
tect and  marine  engineer  in  developing  the  most  economi- 
cal,  efficient  or  the  least  weight  and/or  space  gas  tur- 
bine nuclear  power  plant. 

For  any  ship  propulsion  plant,  it  will  be  neces- 
sary for  the  ship  designer  to  produce  machinery  compon- 
ents which  give  the  best  performance  for  the  least 
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weight,  cost,  space  or  some  other  criteria.   The  subma- 
rine designer  for  instance  would  want  to  optimize  his 
combination  of  machinery  components  in  order  to  achieve 
the  best  performance  for  the  least  weight  and  space. 
The  merchant  ship  designer  on  the  other  hand  would 
direct  his  work  toward  producing  a  plant  that  gives 
the  best  performance  for  the  least  cost.   Optimization 
of  a  propulsion  plant  can  be  performed  for  different 
reasons.   Knowing  the  mission,  size  and  requirements 
of  the  ship,  the  designer  is  able  to  optimize  the  machin- 
ery on  a  basis  that  gives  the  best  power  plant  commen- 
surate with  the  ship's  purpose. 

In  an  optimization  procedure,  the  real  problem 
is  to  reduce  the  number  of  variables  so  that  they  can 
be  conveniently  handled.   The  ideal  situation  is  to 
have  the  weight  or  cost  expression  in  terms  of  one  para- 
meter or  at  most  two  independent  parameters.   Computa- 
tion work  and  time  are  greatly  reduced  with  each  less 
variable.   Also  desirable  is  an  expression  in  closed 
analytic  form  which  will  greatly  simplify  calculations. 
However,  the  functions  to  be  differentiated  will  usually 
not  be  found  in  this  form,  and  the  designer  must  then 
resort  to  a  graphical,  numerical  or  some  other  method 
of  approximate  differentiation.   It  is  important  that 
careful  study  be  given  all  parameters  in  a  weight  or 
cost  function  to  determine  whether  or  not  they  are 
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independent.   If  two  parameters  are  dependent,  the  method 
of  evaluating  the  function  will  be  different  and  gener- 
ally more  difficult  than  if  all  variables  are  independent, 

The  first  step  in  eliminating  the  number  of  vari- 
ables is  to  assign  possible  values  to  those  quantities 
that  are  beyond  the  control  of  the  designer.   Secondly 
the  concept  of  'balance15  is  used  to  reduce  the  number 
of  variables.   The  latter  step  can  be  applied  to  compon- 
ents which  bear  a  specific  relationship  to  each  other. 
For  example,  the  task  of  finding  the  surface  area  of 
two  different  components  for  a  given  power  plant  might 
be  reduced  to  a  single  problem  of  selecting  a  balanced 
combination  of  these  areas.   This  method  would  be  pos- 
sible if  the  effects  of  each  area  on  cost  or  plant 
weight  were  evidenced  by  their  influences  on  the  same 
working  fluid  property.   The  problem  of  selecting  a 
balanced  combination  of  areas  is  done  separate  from  the 
optimization  of  the  plant  as  a  whole.   All  other  compon- 
ents of  the  plant  must  be  sensitive  only  to  this  fluid 
property  and  not  to  the  various  combinations  of  these 
two  areas  which  establish  that  fluid  property.   Once  a 
relationship  between  balanced  components  is  formulated, 
the  units  comprising  the  balanced  group  are  treated  as 
a  single  component  with  a  consequent  elimination  of 
variables. 

In  this  work,  a  number  of  variables  were  eliminated 
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from  consideration  of  step  one.   The  'balance1*  concept 
was  used  to  form  a  relationship  between  the  compressors 
and  the  high  pressure  turbine.   This  step  eliminated 
the  high  pressure  turbine  and  low  pressure  turbine 
pressure  ratios  as  variables.   Thus  once  a  reasonable 
value  for  the  cycle  pressure  loss  is  specified,  the 
remaining  variables  are  flow  rate  (w)  ,  top  pressure 
level  (P) ,  overall  compressor  ratio  (r),  regenerator 
effectiveness  (r\    )  and  temperatures  in  the  cycle. 
Both  cycle  efficiency  and  net  work  per  pound  of  fluid 
flow  are  independent  of  top  pressure  as  shown  in  Sec- 
tion II.   Pressure  is  then  of  importance  only  in  the 
way  it  effects  weight  and  space  of  components.   All 
remaining  temperatures  of  the  cycle  are  dependent  on  r. 
For  a  specific  power  output,  flow  rate  is  a  function  of 
the  pressure  ratio.   With  establishment  of  a  fixed  reac- 
tor inlet  temperature  as  previously  discussed,  regenera- 
tor effectiveness  is  also  a  function  of  pressure  ratio. 
Because  the  reactor  inlet  temperature  is  set,  the  cycle 
pressure  ratio  must  always  be  less  than  the  value  that 
would  make  the  turbine  exhaust  temperature  equal  to  the 
reactor  inlet  temperature.   For  equal  reactor  inlet  and 
turbine  exhaust  temperatures,  a  ti   of  100  per  cent  is 
required.   A  pressure  ratio  giving  a  turbine  exhaust 
temperature  lower  than  reactor  inlet  temperature  would 
require  a  r\      greater  than  100  per  cent  and  hence  presents 


an  impossible  cycle* 

In  view  of  the  above,  there  are  two  independent 
variables  for  a  plant  of  this  type.   All  weight  and  size 
equations  were  developed  in  terms  of  these  two  inde- 
pendent parameters;   pressure  and  pressure  ratio.   Sub- 
sequent solution  of  the  equations  for  a  narrow  range 
of  pressure  ratios  indicated  extreme  variations  in 
total  weight  and  size  of  components. 

The  regenerator  weight  and  size  was  found  to 
dominate  that  of  all  other  components  by  several  orders 
of  magnitude  as  regenerator  effectiveness  approaches 
100  per  cent.   In  order  to  achieve  a  higher  efficiency 
without  having  to  accept  the  large  increase  in  weight 
contributed  by  the  regenerator,  consideration  was  given 
to  the  idea  of  using  an  additional  piece  of  heat  exchang- 
er equipment  between  the  regenerator  and  reactor  inlet, 
A  small  quantity  of  high  temperature  fluid  could  be  re- 
moved from  the  reactor  to  further  heat  the  bulk  of  the 
fluid  from  the  regenerator.   The  extraction  of  this 
fluid  from  the  reactor  might  occur  after  it  has  passed 
up  the  outside  of  the  reflector  in  cooling  the  thermal 
shield  and  vessel  wall  and  before  it  flows  along  the 
surface  of  the  fuel  rod  assemblies.   Another  possible 
extraction  point  for  the  small  quantity  of  high  tempera- 
ture fluid  might  be  in  the  LPT  after  partial  expansion 
or  at  the  HPT  exhaust.   This  bled  fluid  would  then  be 
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recirculated  to  the  additional  heat  exchanger  to  further 
raise  the  temperature  of  the  bulk  of  the  fluid  coming 
from  the  regenerator.   A  considerable  portion  of  the 
regenerator  duty  could  then  be  taken  over  by  the  addi- 
tional heat  exchanger.   The  resulting  lower  regenerator 
effectiveness  would  lower  regenerator  weight,   A  study 
of  such  a  proposal  is  needed  to  indicate  if  the  reduc- 
tion in  regenerator  weight,  caused  by  a  lower  regenera- 
tor effectiveness,  is  more  than  the  weight  increase 
caused  by  the  additional  heat  exchanger  (with  associated 
equipment)  and  the  greater  weights  of  other  components 
(that  might  be  caused  by  the  slightly  larger  flow  rate 
required  in  order  to  keep  the  same  power  output )» 
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II,  PROCEDURE 

A.  General 

The  basic  procedure  followed  was? 

1.  Develop  mathematical  expressions  for  component 
weights  and  sizes  in  terms  of  state  point  properties 
and  system  constants. 

2.  Determine  the  required  system  constants. 

3.  Solve  the  equations  for  a  range  of  operating 
conditions „ 

Steps  1  and  2  are  carried  out  in  this  section 
with  the  detailed  procedures  of  the  development  per- 
formed in  Appendices  A,  B  and  C.   Step  3  is  found  in 
Appendix  D. 

The  plant  components  considered,  consist  of  the 
regenerator,  precooler,  intercooler,  reactor,  accumula- 
tor, piping,  and  rotating  machinery  comprised  of  the 
low  pressure  compressor  (LPC),  high  pressure  compressor 
(HPC),  high  pressure  turbine  (HPT)  and  low  pressure 
turbine  (LPT)a   The  HPT  drives  the  HPC  and  LPC,  and  the 
LPT  provides  shaft  power.   A  low  pressure  turbine  by- 
pass heat  sink  was  not  considered  as  a  plant  component,, 
Preliminary  studies  [13]   indicate  that  direct  bypass- 
ing of  the  low  pressure  turbine  without  an  additional 


■^"Numbers in  brackets  refer  to  references  listed  in 
Appendix  E. 


heat  exchanger  is  feasible.   With  this  system  thermal 
cycling  conditions  appear  to  be  within  the  capabilities 
of  the  different  component  materials.   A  system  requir- 
ing additional  piping,  but  still  no  bypass  heat  sink, 
would  use  a  split-flow  scheme.   A  portion  of  the  by- 
passed fluid  would  be  piped  directly  to  the  precooler 
and  the  remainder  would  proceed  through  the  regenerator 
as  usual.   The  reactor  inlet  temperature  could  be  kept 
more  closely  regulated  with  this  latter  method.   Figure 
11=1  is  a  flow  diagram  of  the  system  indicating  state 
points  and  some  of  the  symbols  used  in  the  development 
of  the  equations.   Figure  II-2  is  the  corresponding 
temperature-entropy  diagram  for  the  cycle. 

The  mathematical  development  that  follows  is  based 
on  the  perfect  gas  relationships  and  on  turbulent  flow 
of  the  working  fluid  in  ducts  and  piping.   In  addition, 
the  velocity  of  fluid  flow  in  the  system  is  considered 
to  be  of  the  order  of  magnitude  where  static  pressure 
and  temperature  are  approximately  equal  to  stagnation 
conditions .   All  gas  properties  are  evaluated  at  the 
mean  fluid  temperature  in  a  particular  component  and 
are  considered  constant  over  the  range  of  temperatures 
encountered  in  that  component 

For  minimum  compressor  work  in  a  two-stage  com- 
pression with  intercooling  to  the  low  temperature  of 
the  cycle  (T-.),  the  pressure  ratio  of  each  compressor 
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Figure  IlZ-1 
Flow  Diagram  for  helium-cooled,  graphite  moderated 
Nuclear  Gas  Turbine  Power  Plant 
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can  be  readily  shown  to  equal  the  square  root  of  the  over- 
all pressure  ratio; 

The  above  condition  also  implies  negligible  pressure 
drop  of  the  fluid  through  the  intercooler. 

Total  compressor  work  can  be  written  for  equal 
compressor  efficiencies  as: 

W.   =  J-J-  (f5  -  1)  (2-2) 

^c 

Total  turbine  work  for  equal  efficiencies  1st 

\     -     VpT6  (1  "  &  <2"3> 

Y  being  the  pressure  loss  factor  in  the  cycle  defined  ass 

T^/TQ   t.   /r.Y   /t>  P2  P-,\v   /Pc  Pc  P0  P,  P-,\v 
Y 6/x8s  __t  _  f  _t  />f_6  _i  JL   _   _6  J2  _2  _I  _2  (?,\ 

"  TWTi  "  T  "  ^J      W  P4  p2>  "  vp5  P4  p8  p9  p2^  C2  4; 

t  is  the  ratio  of  compressor  isentropic  exit  temperature 
to  its  inlet  temperature  and  t,  is  the  ratio  of  turbine 
inlet  temperature  to  its  isentropic  exit  temperature. 
The  factors  in  the  last  term  are  the  fluid  pressure  drops 
contributed  respectively  by  the  reactor,  cold  side  of 
regenerator,  hot  side  of  regenerator,  precooler  and  inter- 
cooler o   The  pressure  drop  being  a  function  of  working 
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fluid s  pressure,  and  component  arrangement  and  geometry- 
is  usually  estimated  by  the  designer  for  each  component . 
The  total  cycle  pressure  drop  and  hence  Y,  for  the  cycle 
or  a  component,  can  be  determined. 

The  net  work  for  the  cycle  is  simply  W,  -  W  . 
Using  (2-2)  and  (2-3)  the  net  work  (W net)  per  pound 
of  fluid  flow  can  be  written  as; 


^  =  "t  if  "  -  fe>  -  f;  <*'5  -  i)       (2-5) 


Heat  into  the  cycle  from  the  reactor  can  be  written 
ass 

Q.      C  (fc   Tc)       Tt  -  T. 

P  !     CP  Tl  Tl 


The  thermal  efficiency  of  the  cycle  is  then 


w. 

I 

'in 


net  (2-7) 


Combining  (2-5),  (2-6)  and  (2-7) 


11 1  T-,  vx   Yt' 


,  -     -  ?  (1 '  k)  -  k  ("5  - 1}    c«: 


Tl 


For  reasons  discussed  in  Section  I,  T^,  T,-,  T1 
r\      and  7]   are  not  generally  subject  to  the  designer's 

Ti         C 
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control o   With  a  reasonable  value  assigned  to  Y,  the  net 
work  and  efficiency  become  functions  of  cycle  pressure 
ratio  only»   By  setting  the  derivative  with  respect  to 
x  equal  to  zero  in  (2=5),  the  optimum  value  of  cycle 
pressure  ratio  (r)  is  found  which  maximizes  cycle  net 
work.   Similarly  the  optimum  value  of  r  can  be  found 
from  (2=8)  which  gives  maximum  cycle  efficiency. 

Bo   Heat  Exchanger  Equipment 

All  heat  exchanger  equipment  considered  is  assumed 
to  behave  as  if  of  the  shell  and  tube  counterflow  con- 
figuration.   In  the  regenerator  the  high  pressure,  cold 
gas  flows  inside  the  tube,  with  the  hotter  gas  flowing 
longitudinally  around  the  tube  interspace.   The  increased 
strength  of  the  tubes  with  internal  pressure  would  neces- 
sitate this  arrangement.   Water  flows  through  the  tubes 
and  gas  around  the  tubes  in  both  coolers.   This  arrange- 
ment would  facilitate  the  required  periodic  cleaning  of 
the  water  side  of  the  coolers.   The  precooler  and  inter- 
cooler  require  large  cooling  water  flows  with  low  water 
velocities  and  low  pump  heads.   Design  factors  for  the 
coolers  follow  marine  steam  condenser  practice, 

I,  Regenerator 

For  a  given  gas,  tube  si*ze  and  tube  arrangement, 
equations  (A-^O)  and  (A-31)  can  be  written  respectively 

-^■Equations  carrying  letter  designations  are  from  the 
corresponding  appendix, 
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as 


1      (T      +  r2    Tu   C,) 


5 


^i  =  %dw(i^)        * 


NdiL  -  %dL  w  (i^; )     ^s  <V  r*W*       <2-10> 

where  C ,  -  —   T~       an^-  "kne  K's  are  constants 


of  proportionality.   Since  T  and  T,  were  defined  as 


Tc  =    2  2  (2~li: 


and   Th  =   8  |  ?  (2=12) 

they  may  be  put  in  terms  of  pressure  ratio  and  known 
cycle  temperatures.   From  the  definition  of  turbine  works 


CP(T6  "  V   =  Vt(T6  -  T8s} 


T8  -  T6[1  ~  V1  "  k>]  (2'15) 

Using  (2=13)  and  (A-4-) ,  (2=11)  and  (2=12)  can  be  written s 

Ts(2^  )  -  7i    T  [1  =  n  (1  -  ^)] 

T  »   -5 x  p  ,-,x   "b. -& LE_-        (2=14) 

c  2  CI  -  tix; 
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xh  "      2(1  -  nx)  ^"li?; 


Substituting  (2-14)  and  (2-15)  into  (2-9)  and 
(2-10),  cross  sectional  and  heat  transfer  area  variations 
become  functions  of  presure  ratio  (r) ,  "n  ,  w,  P  and 

A 

known  quantities;   For  a  specified  maximum  and  mini- 
mum temperature  of  the  working  fluid  (Tfi  and  T,),  a 
fixed  reactor  inlet  temperature  (T,-)  and  an  appropriate 
value  for  turbine  efficiency,  the  area  variations  are 
dependent  on  r,  P,  and  w.  r\     becomes  known  once  the 
pressure  ratio  is  specified,  since  reactor  inlet  tem- 
perature is  fixed. 

Regenerator  weight  may  be  written  as  a  function 
of  the  heat  transfer  area: 

WRQ  =  K-^Nd-L)  +  K2  (2-16) 

where  K,  and  Kp  are  constants  and  (Nd.L)  is  given  by 


(2-10)  with  (2-14)  and  (2=15)  substituted  for  T  and 

v 


2.  Coolers 

For  a  given  gas,  tube  size  and  tube  arrangement, 
equations  (A-32)  and  (A-33)  for  the  precooler  become; 

C\  T  V5     .5 
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T 


NdoL  =  KNdL  w 


1.4 


m 


T 


pc 


8 


where  Z\  T   =  Tq  -  T, 


T 


'pc 


T9  *  Tl 


(2-18) 

(2-19) 
(2-20) 


Us 


ing  equations  (2-13)  and  (A-4) ,  TQ  can  be  expressed  as 


V^x-^Vx  +  Y? 


T, 


1  -  r\. 


(2-21) 


y\  T  ,  the  LMTD  may  be  written  as 


A  T   =  a(z\TA)  =  a 


m 


^AT.  +  AT 

-A  i-  a[T   -  T   ]    (2-22) 
pc   mw 


whe 


re  Z^T.  is  the  arithmetic  mean  temperature  difference 


U-H 


AT. 


and  a  is  dependent  on  the  ratio    i  \    see  Figure  II-3 

v^T 


2  -  distance  along  cooler 
Figure  II- 5 

Cooler  Temoerature-Length  Diagram 
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For  the  precooler  and  intercooler,  a  mean  water  tempera- 
ture (T   )  of  80°F  is  assumed,  based  on  70°F  sea  water 

mw  * 

inlet  temperature  and  a  20°F  water  temperature  rise 

through  each  cooler.   From  [2]^ values  of  V  range  from 

1.00  to  0.710  as   ATi   goes  from  1  to  10. 

-^  e 

With  T  ,  and  T  ~  specified, "a* can  be  approximated  by: 

a  =  1.227  -  .0006  Tg  (2=23) 

for  the  range  of  gas  temperatures  likely  to  be  encoun- 
tered at  the  precooler  inlet, 

With  the  above  simplifications,  (2-19),  (2-20)  and 
(2-22)  may  be  expressed  in  terms  of  pressure  ratio,  r\ 

-A, 

and  known  temperatures  and  efficiencies  by  using  (2-21) 
and  (2-23).   The  resulting  relations  for  (2-19),  (2-20) 
and  (2-22)  may  then  be  substituted  into  (2-17)  and  (2-18) 
giving  rather  long  equations  for  the  cross  sectional  and 
heat  transfer  areas  in  terms  of  flow  rate,  pressure, 
pressure  ratio,  11   and  known  quantities. 

.X. 

For  a  given  working  fluid,  tube  size  and  tube  ar- 
rangement, (A-34)  and  (A-35)  for  the  intercooler  becomes 

Ndo  "  KNd  w(|^]    (*icJ5    4-  <2^ 

NdoL  -  KNdL  *&£)    r*«-V    ^  (2^25) 


[NdL  Wt^Tmj    (Tic)'    p. 8 
=18- 


where    /\  T 


g  =  T2  "  T5  "  T2  '   Tl 


(2-26) 


T. 


T?  +  T? 


T2  +  Tx 


1C 


(2-27) 


T2  represents  the  temperature  out  of  the  LPC  and  can 
be  expressed  from  compressor  work  as: 


T2  =  Tx 


(f5  -  1)  +  1 


Combining  (2-28)  with  (2-26)  and  (2-27): 


AT  =  T. 


-1  r,.5 


T. 


T 


1C 


(t^  -  1) 


_1  (T«5  _  i)  +  2 


The  LMTD  corresponds  to  that  of  the  precooler 


(2-28) 


(2-29) 


(2-50) 


AT  =  a(T.   -  T   ) 

^  m      ic    mw 


(2-3D 


where ;;a^ can  be  approximated  by: 

a  =  1.814  -  .0014  T2  (2-32) 

for  reasonable  values  of  gas  temperature  at  the  inter- 
cooler  inlet.   Substituting  (2-28)  into  (2-32): 


a  =  1.814  -  .0014  T- 


r1  (t-5-  1)  +  1 
c 


(2-33) 
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Combining  equations  (2-30),  (2-31)  and  (2-33): 


/\Tm=  1. 81-. 0014  ^[^-(T^-D+lDj/^C^T-5-!)^]-^^)   (2-34) 


where  5?   was  taken  as  80  +  460  =  540°R. 

mw 

By  placing  (2-29),  (2-30)  and  (2-34)  into  both  (2-24) 
and  (2-25),  the  two  intercooler  area  variations  are 
expressed  in  terms  of  flow  rate,  pressure,  pressure 
ratio  and  known  temperatures  and  compressor  efficiency. 
As  for  the  regenerator,  the  precooler  and  inter- 
cooler weights  may  be  expressed  as  functions  of  the 
heat  transfer  area.   For  the  precooler  the  basic  weight 
equation  is: 

/^T  \1.4 
W   =  K,  wf  rm2     T       /rY8  +  K~         (2-35) 

where  ^\  T  -AT  and  T   can  be  expressed  in  terms  of  P, 
g    m      pc  r 

r,  Tj   and  known  quantities  from  equations  (2-19)  through 
(2-23). 

The  weight  equation  for  the  intercooler: 

_   k     w    fcl&\  /T.  V4  21.     +  v  (2-36) 

ic    "Kl   W    (ATjic     (     ic)  p. 8  2 


W 


where  a  T  ,  /\  T  and  T.   can  be  expressed  in  terms  of  P, 
g     m      ic  ^ 

r  and  known  values  from  equations  (2-29),  (2-30)  and  (2-34). 
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C.  Reactor 

The  reactor  weight  and  size  equations  are  developed 
"by  considering  the  reactor  as  another  piece  of  heat  ex- 
changer apparatus.   The  proposed  MGCR  helium  cooled, 
graphite  moderated  reactor  with  a  heterogeneous  assembly 
of  fuel  rods  [11],  [12],  [14]  is  used  as  the  reference 
reactor. 

For  a  constant  heat  generation  flux  along  the 
coolant  channel  and  radially  over  the  reactor,  the  gas 
temperature  rises  with  a  constant  slope  versus  coolant 
channel  length.   The  coolant  wall  temperature  will  also 
rise  with  the  same  slope  from  (Tc-  +  Z\T)  to  (T^  +  z^T) 

where  AT  is  the  coolant  wall  temperature  minus  the 
temperature  of  the  coolant.   The  maximum  wall  and  cool- 
ant temperature  will  occur  at  the  reactor  outlet  with 
the  outlet  coolant  temperature  in  each  channel  being 
equal  for  uniform  heat  generation  in  the  core.   A  con- 
stant heat  generation  flux  is  substantially  true  for  a 
relatively  small  reactor  core  using  good  reflector  mate- 
rial [1],  [33,  [153  and  [183. 

The  heat  transfer  in  the  reactor  can  be  represented 
by? 

Q  =  H  Ac  AT  (2-57) 

Q  =  w  Cp(T6  -  T5)  (2-37A) 

with  AT  =  maximum  wall  temperature  minus  the  outlet 
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temperature  of  the  coolant.   H  is  given  by  (A-13): 

CT   k'2  (w  C  )*8 

H  "    a  i  « E"2T"  (2~58) 

N.8  D1.8  (pr)-^ 

D  and  N  represent  the  outside  diameter  of  the  fuel  rod 
and  number  of  fuel  rods  respectively.   The  pressure 
drop  through  the  reactor  is  given  by  (A-17).   Tbe  term 

2     I    w  m —  J  in  (A-17)  can  not  be  neglected,  but  is 

equal  to  a  constant  for  uniform  heat  generation  and 
fixed  inlet  and  outlet  coolant  temperatures.   (A-17) 
then  becomes: 

A  jJ£  =  -C2  -|   [4f  §]  (2-39) 

5 

where  for  low  mach  numbers,  static  pressures  and  tempera- 
tures can  be  interchanged  with  stagnation  properties. 
The  pressure  drop  factor  (YR)  for  the  reactor  is: 

XH  .  (^  (2-40) 

Following  the  procedure  used  in  obtaining  equation 
(A-22)  in  Appendix  A,  the  pressure  drop  term  for  the 
reactor  can  be  written: 


yR  =  C$  S=l  m2  (4f  §)  (2-41) 

C,  being  a  constant  determined  by  the  fixed  inlet  and 
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outlet  temperatures  and  the  maximum  wall  temperature  of 
the  reactor.   Using  (A-24)  and  (A-26)  the  Fanning  fric- 
tion factor  and  mach  number  squared  are  represented 
respectively  by: 

/DNC^V2 
t   =   .0438  I   Cp  )  (2-42) 

tf    ,    fkf        1_   (wC    )2        !sE    (l=l     _1_Y  (2-43) 

m  .      m 

T  w  =        6  5      =   constant.      Combining  all   constant   terms ; 


ej 


=  cM   ^ 


PND* 


Substitute  (2-42)  and  (2-44)  into  (2-41): 

^R  "  CM3  ^  K  J  H       p2   N1.8D4.8     (2~45) 

Combine  (2-37)  and  (2-38)  and  write: 

C,  k'2  (wC  )  °8 
Q.        «  i  «   P  ,      (tcNDL)  (2-46) 

AT   =   N'V^rPr)*4 


^T  as  defined  with  equations  (2-37)  is  a  fixed  value 
and  so  (2-46)  can  be  written: 

c  .   <Lf8  (?r)-4  -  (2„47) 
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Multiplying  equation  (2-4-5)  by  (2-47)  and  solving  the 
resulting  expression  for  ND2  : 

(2-37A)  can  be  substituted  for  Q  to  give : 

■*  -  CNB  (¥J"5  wCp  ^       ^-  (2-49) 

The  above  equation  represents  the  cross  sectional  area 
variation  of  the  reactor.   For  fixed  inlet  and  outlet 
temperatures  and  a  given  fuel  element  and  coolant  channel 
configuration,  it  is  dependent  on  the  flow  rate  and  top 
pressure  level  in  the  cycle. 

The  corresponding  heat  transfer  area  (NDL)  equa- 
tion is  obtained  by  solving  (2=47)  for  L,  dividing  by 
D,  making  the  appropriate  substitution  for  (ND2 ) '   and 
multiplying  this  result  with  (2-4-9); 


NDL  -  CNDL 


¥)1,2ffl*V)-6*(iy2cp 


w(p)  (t6-t5)"  £-50) 


In  (2-50),  L  is  the  coolant  channel  length  or  the  core 
height,  H.   Using  the  bare  core  approximation,  the  rela- 
tionship between  height  and  radius  for  a  critical  cylinder 
of  minimum  volume  is  given  by  reference  [4]  as: 

H  =   n2[lo5     R     =   1,85R  (2-5D 
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D  is  fixed  for  a  given  reactor  and  N,  the  number  of  fuel 
rods  (coolant  channels)  is  a  function  of  reactor  cross 
sectional  area:    N  =  CNR2 .   Placing  this  and  (2-51)  into 
(2-50)^ reactor  radius  cubed  becomes: 


R3  =  C 


R 


W,2®-Vr6w8aA 


w 


Reactor  weight  can  be  expressed  as: 
WR  =  ^Rn   R2  H 


^(VV^-52) 


(2-55) 


where  r^  is  the  overall  reactor  density  determined  from 
a  reference  or  similar  reactor.   Using  (2-51)  and  (2-52) 
reactor  weight  becomes: 


WR  "    CE 


—I       {* 


U 


%if2/iVV)-^T  (I 


.2 


w 


78(VT5)1#4     (2~54) 


For  a  given  coolant  and  fuel  element  configuration  it  is: 


wr  =  Ks  ^78  <T6  -  V 


(2-55) 


The  R  and  H  used  above  are  for  the  core  and  do  not 
represent  the  actual  size  of  the  reactor.   The  overall 
reactor  radius  and  height,  exclusive  of  shielding,  can 
be  taken  as  some  constant  multiple  of  the  core  dimen- 
sions for  small  changes  in  core  R  and  H*   From  the  MGCR 
reference  reactor,  the  reactor  vessel  radius  (Ry)  and 
height  (Hv)  can  be  represented  as: 


Rv  =  1.63  (R) 


(2-56) 
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Hy  =  3.51(H)  (2-57) 

D.  Turbo  Machinery 

1.  Compressors 

In  formulating  length  and  weight  equations  for 
the  turbines  and  compressors,  it  was  assumed  that  the 
length  of  each  component  varies  directly  with  the  num- 
ber of  stages  and  that  weight  is  proportional  to  the 
product  of  annulus  area  at  the  low  pressure  end  and 
the  length.   Therefore  by  picking  reasonable  values 
for:   efficiency,  p,   £g  ,  31$  P2  and  £    }   equation 

(B-19)  in  Appendix  B  can  be  written  ass 

NLpc  =  K  Tx    (t"5  -1)  (2-58) 

Since  T, ,  the  minimum  temperature  of  the  cycle  is  fixed, 
the  LPC  length  is: 


L 


'LPC 


=  K£  (t'5  -  1)  (2=59) 


The  annulus  area  given  by  (B-12)  becomes: 

*LPC-Ka   *T  (2-60> 

w  is  the  flow  rate,  r  the  pressure  ratio  and  P  the  top 
pressure  in  the  cycle. 

The  weight  equation,  after  applying  the  condition 
stated  above  to  (2-59)  and  (2-60),  is: 
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WLPC  "  Kw   ¥   ^   "  1)  (2-61) 

Using  (B-ll)  and  (B-15),  the  tip  diameter  at  the  LPC 
low  pressure  end  becomes: 


DIPC  -  (*  +  ()(*§f  •      Kd  (*§f  (2-62> 

The  four  corresponding  equations  for  the  HPC  are 
indicated  below: 

LHPC  =  Kl   (T'5  "  1^  (2"65) 


j-5 

'HPC    "a   P 


V  -   Ka   !Sr  (2-64) 


WHPC  "  Kw  H^      C^5  "  1)  ^2-65) 


The  constants  for  each  of  the  above  equations  are 
not  necessarily  applicable  to  both  the  HPC  and  the  LPC. 
Each  machine  has  a  particular  constant  associated  with 
it. 

2,  Turbines 

As  for  the  compressors,  reasonable  values  of  effi- 

ciences,  0,   _a  ,  0,  ,  0P  and  ±   can  be  assumed  or  ob~ 
U     L  d 

tained  from  the  turbine  manufacturer.   These  factors  will 
change  very  little  over  a  wide  range  of  operating  conditions 
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(pressure,  pressure  ratio,  flow  rate,  etc.)  and  so  the 
size  and  weight  equations  can  be  expressed  similar  to 
those  for  the  compressors. 

For  the  HPT,  (B~$2)  gives  the  number  of  stages 
as: 

Njjpj,  =  K  (t-5  -  1)  (2-67) 

Turbine  length,  annulus  area,  weight  and  tip  diameter 
equations  become  respectively  upon  using  the  basic 
relations  derived  in  Section  III  of  Appendix  B: 

LjjpT  =  Kp  (t-5  -  1)  (2-68) 

AHpT=Ka       |      [1   -  Ht   C(T»5   .   1}]  (2_69) 

[1   _   c(t*5  -   l)]u 

2   Tx 
where   C   = 


T6V< 


WHPT   -  Kw       P * ^2~70) 

[1    -   C(t'5-1)]U 

V  CI   -   C(t'5  -   l)]u      / 


For  the  LPT,    the   equations   are: 

1^  =  V    Ci-i»tcCT^-D3CTT(i-c[T^-i3)-i3  (2^72) 

Yt    [1   -   C(t'5-   1)] 
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ALPT  -  Ka   ¥   [1  "  V1  ~  Y^]  (2"7^ 


WT=K  Is  S H- 1 C. (2-74) 

L  w  P  Yt[1-C(t'5-1)] 


The  constants  associated  with  each  equation  will 
be  different  for  the  HPT  and  LPT  unless  considerable 
similarity  exists  between  the  two  machines. 
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E.  PIPING  AND  ACCUMULATOR 

1.  General 

Both  piping  and  the  accumulator  contribute  signifi- 
cantly toward  the  total  weight  of  the  components  within 
the  containment  vessel.   In  addition,  the  space  occupied 
by  the  accumulator  will  effect  the  size  of  the  contain- 
ment vessel  and  hence  the  amount  of  secondary  shielding. 

Weight  equations  for  piping  and  the  accumulator 
are  developed  in  this  section.   The  constants  associated 
with  these  equations  were  not  determined  due  to  lack  of 
suitable  data  from  a  reference  design.   The  equations 
however  indicate  the  trends  to  be  expected  for  different 
operating  conditions. 

2.  Piping 

Weight  of  piping  is  considered  proportional  to  the 
internal  surface  area  plus  a  constant. 

W  =  K,7tDL  +  K0  (2-76) 

pi       d 

For  negligible  heat  transfer,  the  pressure  drop  term  for 
the  piping  is  expressed  as: 

yp  =   S=l  t   M*  |  (2-77) 

From  (A-24-),  the  Fanning  friction  factor  for  turbulent 
flow  is i 

/DNuC  A  °2 
f  =   •°^(wC-2J  (2~78) 
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Mach  number  squared  is  given  by  (A-26)  as: 


M*       .      C(^l        T_     iK^ll2       1  (2_?9) 

R  p2         K3  D4 

Substitute    (2-78)   and   (2-79)    into    (2-77)    and   solve 

resulting  expression  for  D. 

(uC    )-041  ,„r    /mTV2l/ir  iV62 

v-s-Jm  (wV575(^)  ri1)       (2-80) 

For  a  given  working  fluid,  (2-80)  becomes; 

D  =   K    5L__  (tl)-21  (2_81) 

Since  temperature  in  equation  (2-81)  is  raised  to  only 
the  0.21  power,  temperature  does  not  effect  the  diameter 
(weight)  significantly.   As  a  preliminary  conservative 
weight  estimate,  the  maximum  temperature  in  the  cycle  is 
used  with  (2-81).   For  greater  accuracy,  the  piping  length 
could  be  divided  into  two  halves,  a  hot  and  a  cold  section. 
Equation  (2-81)  can  now  be  written  as: 

w'375     P1 
D  =  K   *-_-  (L)-21  (2-82) 

-P° 
Combining  (2-76)  and  (2-82): 

W   =   K  w°41   L  lo21  +  K2  (2-8$) 

Considering  the  length  of  piping  within  the  con- 
tainment vessel  as  a  fixed  quantity  for  similar  ships, 

(2-8$)  becomes: 
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w-375 


The  constants  in  (2-84)  were  not  determined  since 
suitable  information  from  a  reference  design  concerning 
length  and  weight  of  piping  was  not  available. 

3.  Accumulator 

The  accumulator  weight  is  expressed  as  a  function 
of  the  total  weight  of  the  working  fluid  in  the  cycle 
and  the  allowable  stress  in  the  wall  material.   Consider- 
ing a  long  cylindrical  container,  the  surface  area,  upon 
neglecting  the  end  plates,  is; 

A   =  2uRL  (2-85) 

The  volume  of  the  working  fluid  is  given  by: 

V  =   tcR2L  (2-86) 

V  =  m  S£  (2-87) 

P 

where  m  represents  the  total  mass  of  the  working  fluid 
within  the  cycle.   Weight  of  the  accumulator  is  simply 
the  surface  area  times  the  product  of  wall  thickness 
and  material  density . 

WA  =   27iRLt/^  (2-88) 

Using  simple  hoop  stress  analogy,  wall  thickness  becomes: 

t   .   (P~Fa  )  S  (2-89) 
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where   Pa  is  the  external  (atmospheric)  pressure  and  P 
is  the  pressure  of  the  fluid  within  the  accumulator. 
For  rapid  changes  of  load,  the  pressure  of  the  fluid 
within  the  accumulator  must  be  greater  than  the  top 
pressure  of  the  cycle.   Equation  (2-87)  indicates  a 
high  pressure  is  also  desirable  from  the  standpoint 
of  reducing  the  size  of  the  container.   To  facilitate 
rapid  changes  in  plant  load  [10] ,  the  accumulator  is 
arranged  in  two  parts:   a  high  and  low  pressure  sec- 
tion.  When  reduction  in  loads  occur,  gas  is  bled  from 
the  HPC  discharge  to  a  low  pressure  accumulator  or 
receiver.   As  the  low  pressure  flask  fills,  it  is 
emptied  by  a  piston  type  compressor  which  compresses 
the  gas  to  a  pressure  considerably  above  the  top  press- 
ure in  the  cycle.  The  small  piston  compressor  discharges 
to  the  HP  accumulator.   As  loads  again  increase,  the  HP 
accumulator  discharges  into  the  high  pressure  line 
between  the  HPG  and  the  regenerator.   Equating  (2-86) 
and  (2-87)  and  solving  for  the  radius  squared  (R2): 

R*    =    Bgg  (2-90) 

Combining    (2-88),    (2-89)    and   (2-90),    the   weight   of 
accumulator  becomes: 

.        2^mRT        £faj  (  , 

A  Op 
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Since  the  temperature  of  the  gas  to  be  put  in  the  ac- 
cumulator will  be  of  the  order  of  200  -  300°F  (the  HPC 
discharge  gas  temperature),  the  accumulator  material 
need  only  be  a  carbon  steel.   Therefore  for  a  given 
gas  and  a  moderate  storage  temperature,  equation 
(2-91)  reduces  to: 


WA   =   K  m  (F~Pa)  (2-92) 


Because  the  pressure  within  the  accumulator  will  be 


several  orders  of  magnitude  larger  than  atmospheric 
-] 

P 


P-P 

pressure, &.       ^^1  and  (2-92)  becomes s 


WA   =   K   |  (2-93) 

Maintaining  a  constant  hoop  stress  for  a  fixed  cylinder 
radius  means  that  the  pressure  to  wall  thickness  ratio 
must  be  a  constant.   With  these  restrictions,  accumulator 
weight  is  written: 

WA        KA  m  (2-94) 
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F,  Evaluation  of  Constants 

The  mathematical  development  of  the  weight  and 
space  equations  necessitated  the  introduction  of  a  num- 
ber of  unknown  constants.   Some  constants  depend  upon 
materials  used,  some  are  proportionality  constants  and 
others  are  constants  of  design. 

To  solve  the  equations,  it  is  necessary  to  evalu- 
ate each  constant.   Accurate  evaluation  of  the  constants 
requires  a  thorough  investigation  of  similar  plants, 
with  note  taken  of  the  specifications,  weight  state- 
ments, heat  balances  and  other  information. 

The  equations  developed  under  Section  II  and 
Appendices  A  and  B  are  general  in  nature  and  not  re- 
stricted to  any  particular  working  fluid  or  type  of 
ship.   Only  in  the  evaluation  and  application  of  the 
constants  does  the  problem  become  one  of  a  specific 
type  or  design. 

Since  no  shipboard  nuclear  closed  cycle  gas  tur- 
bine plant  exists,  available  data  for  evaluating  the 
constants  was  limited.   The  author  relied  mainly  on 
information  from  the  proposed  MGCR  propulsion  plant 
presented  in  references  [11],  [12],  [14],  [151  and, 
where  applicable  on  information  from  existing  open 
cycle  gas  turbine  plants.   The  heat  balance,  specifica- 
tions, and  weight  estimates  for  the  helium-cooled, 
graphite  moderated  MGCR  plant  allowed  satisfactory 
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evaluation  of  the  constants  for  all  components  except 
piping  and  the  accumulator. 

The  evaluation  of  these  constants  and  sample 
calculations  are  given  respectively  in  Appendices  C 
and  D. 
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Ill,  RESULTS 

A.  Cycle  Conditions 

Original  results  were  obtained  in  the  form  of  data 
from  hand  calculations  of  the  developed  equations  with 
reference  to  the  helium  cooled,  graphite  moderated  MGCR 
power  plant.   Calculations  were  made  by  varying  the 
compressor  pressure  ratio  (r)  for  two  different  reactor 
inlet  temperatures  (T,-):   T,-  =  700°F  and  T,-  =  791°F. 

Table  III-A  lists  the  fixed  values  used  for  both  sets 
of  calculations. 

Table  III-A 
Fixed  Cycle  Conditions 


LPT  (net  power)  output   =21,500HP   Pressure  drop=— 


Working  Fluid  =Helium 

Top p ressure  (P) 

Maximum  Cycle  Temperature 

(T6)  =1300°F 

Minimum  Cycle  Temperature 


=  .86 


Coolant  Water  Inlet 
Temperature (Twl)    =70°F 

=1000  psia  Coolant  Water  Outlet 

Temperature (T  2)    =90  F 


(T2) 


T^   =  100"F 


Turbine  Efficiency 

(V 


=  .90 


Compressor  Efficiency 
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B.  Graphs 

The  data  is  presented  in  graphical  form  in  accordance 
with  Table  III-B. 

Table  III-B 
Graphs 

Plot  Figure 

Total  Weight  and  Efficiency  vs  Pressure  Ratio  III-l 

Regenerator  Weight  vs  Pressure  Ratio  III-2 

Precooler  and  Intercooler  Weight  vs  Pressure  Ratio   III-3 

Reactor  Weight  vs  Pressure  Ratio  III-4 

Rotating  Machinery  Weight  vs  Pressure  Ratio  III-5 

Regenerator  Dimensions  vs  Pressure  Ratio  III-6 

Precooler  Dimensions  vs  Pressure  Ratio  III-7 

Intercooler  Dimensions  vs  Pressure  Ratio  III-8 

Reactor  Dimensions  vs  Pressure  Ratio  111=9 

Rotating  Machinery  Length  and  Low  Pressure 

Turbine  Exhaust  Diameter  vs  Pressure  Ratio  111-10 

Flow  Rate  (lb/HP-HR)  vs  Pressure  Ratio  III-ll 
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Figure   III-2 


UJ 

cr 

5s 

k 

5 

<x 

Uj 

a 

K 

i 

3: 

U] 

4 

h 

UJ 

s 

h 

UJ 

-J 

X 

=5 

0 

""* 

Q- 

UJ 
UJ 

o 
K 

6 

UJ 

UJ 

c 

\J0 


o 
10 


o 

Cir 


s 


o 
o 


O 


o 


5/V  O-L 


A0- 


Figure   I I I- 3 


DC       * 
uj 


UJ 


UJ 

or 
I- 

d 

or 
uj 
a. 

%. 

uj 
h 

UJ 

a: 
o 

h 
o 

Ul 

Q: 


h^ 


o 


cr 

UJ 

a: 

v\ 
UJ 

or 
a. 


o 
or 


>o 


o 
^1 


S  A/O 


_G 


=41- 


Figure    I I 1-4 


5 


0 
h 

U 

uj 
K 


uJ 

ec 
rs 
K 
Q 
it 
UJ 
Q. 

UJ 
UJ 

<; 

Q: 

o 

K 
U 
«t 

UJ 

ft 


K^ 


o 


vO 


-or 

o 


o 

o 


O 


S/VOJ. 

A2- 


Figure   II 1-5 


o 
<5h 


It 

UJ 

VJ 

K 
O 
t 


5 


UJ 

o: 

i/l 
1/1 
UJ 

o: 
a. 


o 

Or 


O 


X 


w? 


uO 


■0 


s/vo 


43- 


o 
no 


53 


Figure    III 

LENGTH 

o 

~6 
(ft) 

O 

o 

CS 


•9 


0 

Ul 


kfl 


Q         S 
QC    Q  UJ 


in 

0 
<X 
U 


K 

Ct 
UJ 

01 


h 

cr 
o 

h- 
u 

ct 


hu 


in 


kO 


o 


o 


h 

or 


id 
or 

en 
t/) 
UJ 
Ct 

a 


o 


<3~ 


T 


s 


CO 


o 

OS 


-44- 


-O 
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IV.    DISCUSSION  OF  RESULTS 

A.  General 

The  results  obtained  in  this  thesis  are  based  on 
the  Maritime  Gas-Cooled  Reactor  (MGCR)  plant  and  are  ap- 
plicable only  to  a  plant  of  this  type.   Accuracy  of  re- 
sults are  dependent  on  the  accuracy  with  which  system 
constants  were  evaluated. 

Hand  calculations  were  performed  to  solve  the 
developed  equations  for  a  range  of  operating  conditions. 
Results  are  presented  graphically  in  Section  III. 

For  all  components,  the  sharp  increase  in  weight 
or  cross  sectional  area  indicated  by  the  curves  at  low 
pressure  ratios  is  a  result  of  the  large  flow  rates  re- 
quired at  these  ratios.   The  pressure  ratios  which  give 
a  mathematical  regenerator  effectiveness  of  100  per  cent 
are  approximately  3»1  and  3«8  respectively  for  the  high 
and  low  reactor  inlet  temperatures.   A  higher  pressure 
ratio  than  3.1  for  the  cycle  with  T^  =  791°F  would  pre- 
sent an  impossible  plant  since  a  regenerator  effective- 
ness greater  than  100  yb  is  required.   A  similar  state- 
ment applies  to  the  cycle  with  the  lower  reactor  inlet 
temperature. 

B.  Total  Weight 

Total  weight  and  efficiency  are  plotted  against 
pressure  ratio  for  two  different  reactor  inlet 
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temperatures  in  Figure  III-l.   Minimum  weight  occurs  at 
a  pressure  ratio  of  about  2.5  and  gives  an  efficiency- 
equal  to  30.5  °X>  for  TV  =  791  °F.   For  a  reactor  inlet 
temperature  of  700°F  a  slightly  higher  minimum  weight 
occurs  at  a  pressure  ratio  equal  to  3.15  with  efficiency 
again  equal  to  30o5  /o.   For  pressure  ratios  greater 
than  that  for  optimum  weight,  the  weight  rises  sharply 
to  infinity  due  to  the  regenerator  weight.   This  is  be- 
cause of  the  large  increase  in  regenerator  heat  transfer 
area  required  as  turbine  exhaust  temperature  approaches 
reactor  inlet  temperature.   The  MGCR  proposed  plant  is 
designed  for  a  pressure  ratio  of  2.6,  a  value  slightly 
above  that  for  optimum  weight,  but  giving  an  efficiency 
of  32  %. 

Figure  III-l  indicates  that  a  greater  gain  in  ef- 
ficiency is  possible  with  small  changes  in  pressure  ratio 
above  that  ratio  which  gives  optimum  weight.   Although 
lowering  the  reactor  inlet  temperature  gives  a  flatter 
optimum  weight,  overall  efficiency  is  significantly  re- 
duced since  more  heat  must  be  put  into  the  cycle  for  the 
same  reactor  outlet  temperature.   In  addition,  the  opti- 
mum weight  for  T[-  =  700°F  is  larger  than  that  for  a 
higher  inlet  temperature  because  of  the  increased  reac- 
tor weight.   This  greater  weight  is  caused  by  the  increase 
in  temperature  rise  associated  with  a  low  inlet  temperature 
as  indicated  by  equation  (2-55) •   The  lower  inlet 
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temperature  would  also  tend  to  add  to  the  complexity  of 
reactor  construction. 

C.  Regenerator 

Regenerator  weight  reaches  an  optimum  at  a  lower 
pressure  ratio  than  that  for  total  weight,  and  rapidly 
approaches  infinity  as  T|   goes  to  100  yb.  From  Figures 

•A. 

III-l  and  III-2,  it  is  noted  that  the  regenerator  is  the 
controlling  element  in  total  weight  at  high  pressure 
ratios.   The  designer  desires  a  high  pressure  ratio  in 
order  to  achieve  more  efficiency,  but  without  the  large 
increase  in  weight  that  must  be  accepted  as  indicated 
by  these  figures. 

A  100  yo  regenerator  effectiveness  can  never  be 
realized,  and  in  reality  the  designer  must  be  willing 
to  pay  heavily  in  terms  of  weight  and  space  for  a  regen- 
erator effectiveness  as  low  as  90  yo.  When  one  con- 
siders the  weight  and  efficiency  curve  for  T,-  =  791  F, 
it  is  worthy  to  note  that  efficiency  can  be  increased 
5  points  by  going  from  a  pressure  ratio  of  2.5  to  3»0. 
However,  for  this  increase  in  efficiency,  total  weight 
is  more  than  doubled.   Undoubtedly  this  sizeable  gain 
in  efficiency  could  not  be  realized  for  such  a  weight 
penalty. 

A  possible  solution  in  obtaining  part  of  this  in- 
crease in  efficiency  with  a  less  drastic  increase  in 
weight  might  be  to  use  an  additional  piece  of  heat 
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exchanger  equipment  between  the  regenerator  and  reactor 
inlet.   A  small  quantity  of  high  temperature  working 
fluid  could  be  taken  from  the  reactor.   This  extraction 
might  occur  after  the  fluid  has  passed  up  the  outside 
of  the  reflector  in  cooling  the  thermal  shield  and  ves- 
sel wall  and  before  it  flows  along  the  surface  of  the 
fuel  rod  assemblies.   Another  extraction  point  for  the 
small  quantity  of  high  temperature  fluid  might  be  in 
the  LPT  after  partial  expansion  or  at  the  exhaust  of 
the  HPT.   This  extracted  or  bled  fluid  would  then  be 
recirculated  to  the  additional  heat  exchanger  to  further 
raise  the  temperature  of  the  bulk  of  the  fluid  coming 
from  the  regenerator.   A  considerable  portion  of  the 
original  regenerator  duty  would  then  be  taken  over  by 
the  additional  heat  exchanger  (pre-heater) .   The  result- 
ing lower  regenerator  effectiveness  would  lower  regenera- 
tor weight.   The  heat  transmission  in  this  proposed  pre- 
heater  would  be  very  high  since  both  sides  of  the  addi- 
tional heat  exchanger  will  be  near  the  top  pressure  in 
the  cycle.   With  the  high  rate  of  heat  transmission  pos- 
sible under  these  conditions,  it  may  be  possible  to  keep 

the  size  of  the  pre-heater  quite  small.   Although  this 

ft  w 

small  quantity  of  bled  fluid  would  again  enter  the 

reactor  at  a  higher  temperature  than  the  bulk  of  the 
gas,  the  thermal  inertia  of  the  reactor  and  the  large 
quantity  of  fluid  should  prevent  any  harmful  cyclic 
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thermal  stresses.   In  order  to  produce  the  same  power, 
a  larger  flow  rate  of  fluid  is  required  with  the  addi- 
tional heat  exchanger.   In  addition,  a  small  pump 
(compressor)  would  be  required  to  raise  the  pressure 
of  the  bled  gas  to  reactor  inlet  pressure.   A  study 
of  such  a  proposal  would  indicate  if  the  reduction  in 
regenerator  weight,  caused  by  a  lower  regenerator  effec- 
tiveness, is  more  than  the  weight  increase  caused  by  the 
additional  heat  exchanger  (with  associated  equipment) 
and  the  greater  weight  of  other  components  associated 
with  the  slightly  larger  flow  rate.   Plant  complexity 
with  this  system  may  be  such  that  would  not  warrant  its 
use. 

Figure  III-6  indicates  the  variation  of  regenera- 
tor length  from  zero  to  infinity  as  pressure  ratio  in- 
creases and  more  heat  transfer  area  is  required.   The 
curves  indicating  length  in  this  figure  are  for  each 
unit  of  the  3  unit  regenerator  component. 

D.  Precooler 

Figure  III-3  shows  for  a  pressure  ratio  of  2  and 
larger  that  precooler  weight  has  a  slight  negative 
slope.  This  is  to  be  expected  since  there  is  less  and 
less  heat  transfer  area  required  as  the  turbine  exhaust 
temperature  is  lowered  with  increasing  pressure  ratio. 
For  the  lower  reactor  inlet  temperature,  the  precooler 
weight  increases.   This  effect  is  due  to  the  increased 
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heat  transfer  area  required  of  the  precooler  as  regenera- 
tor effectiveness  decreases. 

Precooler  length  rises  with  a  decreasing  pressure 
ratio  as  shown  in  Figure  III-7.   This  would  be  neces- 
sary since  there  is  little  temperature  drop  across  the 
turbines  and  the  fluid  must  be  essentially  cooled  from 
the  reactor  outlet  temperature  to  the  fixed  low  tem- 
perature of  the  cycle.   A  significant  increase  in 
length  is  required  for  T^  =  700°F  since  the  precooler 
must  perform  more  heat  transfer  as  11  decreases. 

.A. 

E.  Intercooler 

The  intercooler  weight,  as  shown  by  Figure  III-3, 
increases  slightly  at  higher  pressure  ratios.  As  the 
LPC  discharge  temperature  is  raised  with  pressure  ratio, 
more  heat  transfer  area  is  required  to  reduce  this  tem- 
perature to  the  fixed  low  temperature  of  the  cycle. 
There  is  no  change  in  intercooler  heat  transfer  area 
(weight)  with  variation  in  reactor  inlet  temperature. 
The  intercooler  must  always  cool  the  gas  from  the  LPC 
discharge  temperature  to  the  cycle's  bottom  temperature. 

Figure  III-8  indicates  that  intercooler  length  in- 
creases with  pressure  ratio.   For  low  ratios,  little 
heat  transfer  area  is  required  to  bring  the  gas  back  to 
the  cycle's  minimum  fluid  temperature. 

F.  Reactor 

Reactor  weight  increases  sharply  at  low  pressure 
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ratios  since  a  bulky  heat  generating  unit  is  required 
to  raise  the  temperature  of  the  large  flow  rates. 
Since  reactor  weight  is  directly  proportional  to  flow 
rate,  weight  will  decrease  with  less  flow.   It  is  also 
worthy  of  note  that  for  a  drop  of  91°F  in  reactor  in- 
let temperature,  a  weight  increase  of  approximately 
35  tons  results  at  a  pressure  ratio  of  3»0-   This  is 
to  be  expected  since  reactor  weight  as  given  by  (2-55) 
is  proportional  to  temperature  rise  to  the  1.4  power. 
Reactor  vessel  diameter  and  height  curves  have  the  same 
form  as  the  reactor  weight  curves.   The  lower  inlet  tem- 
perature requires  larger  reactor  dimensions  since  equa- 
tion (2-52)  shows  that  radius  cubed  is  proportional  to 
the  temperature  rise  in  the  reactor  to  the  1.4  power. 

G.  Rotating  Machinery 

As  indicated  by  Figure  III-5,  rotating  machinery 
has  a  minimum  weight  for  a  pressure  ratio  of  slightly 
less  than  2.0.  At  low  pressure  ratios,  the  large  an- 
nulus  area  required  causes  an  increase  in  the  weight. 
At  higher  pressure  ratios,  the  weight  increase  is  not 
as  sharp  since  length  of  turbo  machinery  has  less  ef- 
fect on  weight  than  the  diameter  or  annulus  area. 

Rotating  machinery  length  being  a  strong  func- 
tion of  the  number  of  stages  also  increases  sharply 
with  pressure  ratio.   Few  stages  are  required  at  low 
ratios  and  many  stages  are  needed  at  higher  pressure 
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ratios. 

Tip  diameter  for  only  the  low  pressure  turbine 
exhaust  is  indicated  in  Figure  111-10  since  its  dimen- 
sions will  be  larger  than  that  of  the  other  5  rotating 
machines. 

Although  the  rotating  machinery  weight  is  sig- 
nificant, the  length  of  the  unit  is  likely  to  have  a 
greater  influence  on  design  of  the  plant.   A  machine 
of  great  length  would  necessitate  a  larger  contain- 
ment vessel  with  a  corresponding  increase  in  shielding 
weight.   An  in-line  arrangement  of  compressors  and  tur- 
bines would  probably  result  in  lower  total  rotating 
machinery  weight  and  fewer  shaft  sealing,  bearing  and 
support  problems  than  an  arrangement  where  the  low  pres- 
sure turbine  is  in  a  separate  casing.   The  latter  ar- 
rangement would  permit  a  shorter  overall  length  by 
placing  the  LPT  beside  the  HPT  and  compressor  unit. 
This  arrangement  would  have  a  slight  detrimental  effect 
on  efficiency  because  of  losses  in  piping  between  the 
HPT  and  LPT,  but  may  be  necessary  if  containment  space 
is  limited. 
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V.  CONCLUSIONS 

An  analysis  of  the  results  supports  the  following 
conclusions: 

1.  Optimum  weight  for  the  reference  plant  occurs  at  a 
value  of  pressure  ratio  that  is  lower  than  the  value 
that  gives  best  efficiency. 

Efficiency  at  optimum  weight  (pressure  ratio  of 
2.5)  is  30*5  /o .  This  compares  with  a  maximum  possible 
efficiency  of  35.5  ^o  at  a  pressure  ratio  of  3.1 
(corresponding  to  a  regenerator  effectiveness  of  100  /q) 
for  a  reactor  inlet  temperature  of  791°F. 

For  a  reactor  inlet  temperature  of  700°F,  the  ef- 
ficiency at  optimum  weight  (pressure  ratio  of  3*15)  is 
30»5  /o.  The  maximum  possible  efficiency  is  33»3  /o 
occurring  at  a  pressure  ratio  of  approximately  3»8 
(corresponding  to  a  regenerator  effectiveness  of  100  /o). 

2.  With  a  plant  of  this  type,  the  designer  must  be  will- 
ing to  accept  a  heavy  weight  penalty  for  slight  effici- 
ency gains  as  pressure  ratio  is  increased  above  the 
optimum  weight  value  to  the  ratio  that  corresponds  to  a 
regenerator  effectiveness  of  100  yo.  Table  V-A  indi- 
cates the  changes  in  efficiency  and  total  weight  for 
small  increments  of  pressure  ratio  above  that  for 
optimum  weight. 


•58- 


Table  V-A 
Changes  in  Efficiency  and  Total  Weight  vs  Pressure  Ratio 

Pressure  ratio  increment  ^n(°7o)  /\  W(tons) 

Reactor  inlet  temperature  =  791°F 

2.5  -  2.6  1.0  5 

2.6  -  2.7  1.0  7 

2.7  -  2.8  0.8  15 

2.8  -  2,9  0.8  55 

2.9  -  3.0  0.9  290 

Reactor  inlet  temperature  =  700  F 

5.15  -  3.25  .5  4- 

3.25  -  3-35  .4-  5 

3.35  -  3.^-5  .4-  10 

3.4-5  -  3.55  .5  16 

3-55  -  3.65  .4-  4-0 

3.65  -  3.75  .5  165 

3.  A  lower  reactor  inlet  temperature  gives  a  flatter  but 
higher  optimum  weight.   A  significant  rise  in  reactor 
weight  results  from  a  lower  reactor  inlet  temperature. 

4.  The  regenerator  is  the  controlling  component  in 

this  cycle.   As  heat  transfer  requirements  go  up  (regen- 
erator effectiveness  approaches  100  /b),  regenerator 
weight  rises  sharply  and  limits  the  gain  in  efficiency 
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that  is  otherwise  possible  with  a  larger  pressure  ratio, 

5.  The  precooler  and  intercooler  weights  and  dimensions 
experience  a  minimum  over  a  wide  range  of  possible 
operating  pressure  ratios  and  do  not  strongly  influence 
the  total  weight. 

6.  Length  restrictions  imposed  by  available  space  are 
much  more  likely  to  influence  rotating  machinery  design 
than  the  weight  of  this  component. 
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VI.  RECOMMENDATIONS 

The  following  recommendations  are  made: 

1.  Since  most  of  the  developed  equations  show  that  weight 
is  inversely  proportional  to  pressure  (to  a  power),  the 
tendency  would  be  to  increase  pressure  level  to  reduce 
weight  and  plant  size.   However,  maximum  pressure  limi- 
tations would  be  set  from  the  high  temperature,  pressure 
effects  on  materials.   For  a  more  complete  analysis  an 
investigation  should  be  made  to  determine  the  effect  of 
pressure  on  the  reactor,  piping  and  regenerator  weight. 

2.  Equations  should  be  refined  to  give  greater  accuracy. 
Accuracy  can  be  improved  as  plants  of  this  type  are  de- 
signed and  built.   A  survey  of  these  plants  for  compon- 
ent weight  and  size  data  will  contribute  to  a  more  pre- 
cise evaluation  of  system  constants.   Such  a  survey 
would  provide  information  for  proper  evaluation  of  the 
piping  and  accumulator  constants. 

3.  Shielding  accounts  for  a  major  portion  of  total 
plant  weight,  with  the  secondary  shielding  weight  being 
largely  a  function  of  size  and  arrangement  of  components 
within  the  containment  vessel.   A  thorough  study  of  size 
and  arrangement  of  components  is  necessary  to  ensure  a 
minimum  shielding  weight.   Reference  [93  gives  a  suit- 
able method  for  the  preliminary  evaluation  of  shielding 

weight. 
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4.  The  regenerator  weight  and  size  dominates  that  of  all 
other  components  by  several  orders  of  magnitude  as  regen- 
erator effectiveness  approaches  100  /o.  An  investigation 
as  to  the  feasibility  of  recirculating  a  small  quantity 
of  fluid  from  the  reactor  after  it  has  cooled  the  thermal 
shield  or  from  the  turbine  (after  partial  expansion)  through 
an  additional  heat  exchanger  (pre-heater) ,  located  between 
the  regenerator  and  reactor  inlet,  is  necessary.   Such  a 
study  would  indicate  if  the  reduction  in  regenerator  weight, 
caused  by  a  lower  regenerator  effectiveness,  is  more  than 
the  weight  increase  caused  by  the  additional  heat  exchanger 
and  the  greater  weight  of  other  components  associated  with 
the  slightly  larger  flow  rates. 
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VII.    APPENDIX 
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APPENDIX  A 
HEAT  EXCHANGER  EQUATIONS  DEVELOPMENT 

I.  General 

Development  of  equations  for  the  heat  exchanger 
cross  sectional  and  heat  transfer  area  variation  are 
based  on  reference  [17]. 

Pressure  drop  of  the  gas  due  to  entry  and  exit 
effect  in  the  exchangers  are  considered  negligible. 

II.  Regenerator 


Z  -  distance  along  regenerator 
Figure  A-l 

Regenerator  Temperature-Length  Diagram 

The  rate  of  heat  transfer  from  the  hot  fluid  to 
the  cold  fluid  for  equal  flow  rates  and  heat  capacities 
on  each  side  may  be  written  in  the  forme 

(1-1) 


Q  =  wCp   (T5  -  T4) 

Q  =  %   (T8  -  Tg) 


Q  =   UA  AT 
^      s   m 


(1-2) 
(1-3) 
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where  ATm  =  js  T±      =  aT£  is  the  log  mean  temperature  dif- 
ference across  the  counter  flow  regenerator  shown  in  Figure 
A-l. 

Regenerator  effectiveness  is  defined  as: 
m      T   -  T    _   T8  -  Tq 

Combining  equations  (A-l),  (A-2),  (A-3)  and  (A-4) : 

AsU=  Af-=wCp/^   )  (A-5) 

m      y  I  1  -ti  I 

The  overall  coefficient  of  heat  transfer  is  written  as 


JL    +£nUo/&±)    +    i  (A_6) 


UA„    H,  A  .    -t-t- 

s     h  sh    k  ircL      H  A 

c  sc 

For  a  gas-to-gas  heat  exchanger,  the  thermal  resis- 
tance of  the  tube  wall  may  be  neglected  in  comparison  with 

the  film  resistances  and  (A-6)  may  be  written: 
,     ,   A      .      .,     /    H    d.  \ 

*■  k^+  k  m  k    i1  +  \    t)       (A"7) 

h  sh     c      c   \     h    o  ' 

where  A  =  A   ,  the  heat  transfer  area  inside  the  tube. 
s    sc ' 

Reference  [5]  expresses  the  heat  transfer  coefficient 

of  a  turbulent  fluid  as: 

H  =   .023   s*  (Re)°8  (Pr)'4  (A-8) 

e 
The  hydraulic  diameter  outside  of  the  tubes  (hot 

side)  becomes 

deh=   dQq  (A-9) 

where  q  for  triangular  tube  spacing  is: 

V   O 

The  equivalent  outside  flow  diameter  per  tube  can  be 

written  as: 

dx  =   cL0(q)'5  (A-10) 
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For  the  cold  side,  Reynolds  number  can  be  written; 

Re)c  -   ^wCD  (A-ll) 

*d  i  H*cCp 

and  for  the  hot  side : 

4w  C 

Re>h  ■    srfe  <A-12> 

Q  Hi  p 

Substituting  (A-ll)  into  (A-8)  and  noting  that 

d   =  d.  for  the  cold  side: 

nOQ  k  *  (wC  )°8 

Hb    ■     1  Cl  8  f  \     4  (A"1?) 

1         c 

(A-8)  becomes  for  the  hot  side: 

2  /ti.„   x.8 


Hv 


028  V   <WV  (A-14) 


r8  d^'8  q(PrJh^ 
Combining  equations  (A-7),  (A-13)  and  (A-14-): 

.028  (WC  ),8k*2  i   / 

p    c   *  ' 

/?    kc  °2  (Pr)  °4 
where  „£  =   ,—    v   yh    *££      1  for  the  temperature  range 

consideredc   The  overall  heat  transfer  coefficient  then 

becomes: 

V  =   •028(WCp)-8ke-'     /  \  (A_16) 

N.8d .1.8(pr)  .4   li  +  qfl°]-8/ 

1       c  \a. f 

1 

The  pressure  drop  of  a  gas  flowing  through  a  tube 
with  combined  friction  and  heat  transfer  can  be  expressed 
for  low  mash  numbers  [7]  as: 
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dP      KM2  /   L 


O    ' 


where  f  is  the  Fanning  friction  factor  defined  as 

shear  stress     T   .,  ,      . , 
a .   In  the  regenerator  and  coolers,  the 

v    &o   /t   -  T  \ 
term    i  I  S    °  )  is  negligible  compared  to  1  and  is 
2\         TQ  J 

dropped.   With  the  lowmach  number  likely  to  be  encountered 
in  the  ducting  and  heat  exchanger  equipment,  static  pres- 
sures and  temperatures  can  be  substituted  for  stagnation 
pressures  and  temperatures  with  little  loss  of  accuracy. 
With  these  simplifications,  (A-17)  can  be  integrated  to 
give  for  the  cold  side: 

P,        M  2        T 
i  n  «*  -   -K  -f  4  f   §-  (A-18) 

and  for  the  hot  side: 


P 

'8  "  ~o 

From  equation  (2-4),  the  pressure  drop  factor 
(YRG)  contributed  by  the  regenerator  is: 

*HG   "  (^J  <A"2°> 

Take  antilog  of  (A-18)  and  (A-19)  and  combine  with 

(A-20)  to  obtain: 

/   ^    /M  2f    Mj  f  ^\\t 

(A-21) 


Y 


exP  -  2  4Ll"a^ + 

-   v       i 


d  q 
0H 


YRG  is  a  constant  for  the  working  fluid  used  and  for 
a  particular  top  pressure  in  the  cycle. 

After  taking  the  logarithm  of  (A-21): 
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y™  is  a  constant  since  YRG,  as  noted  above,  is  constant. 

The  relationship  for  the  Fanning  friction  factor 
(f)  for  turbulent  flow  can  be  given  as: 

f  -  *MS  2  (A-25) 

(Re)'2 

Substituting  (A-ll)  and  (A-12)  into  (A-23): 

.2 

(A-24) 


(A-25) 


From  the  mach  number  definition  and  the  continuity 
relation,  M2  for  the  cold  and  hot  side  may  be  written 
respectively  as: 

*■  (?y  k(^y  i  (^  *)2   u-26) 

*-  (?)'  fc«V  i(¥  Wo)'     (A"27) 


fc  - 

.04* 

\  wcp 

fh  = 

.044 

I WCP 

where  T   =    T5  +  T4   and  T.     =   T8  +  ^)   • 
"2 h    2 


c 


Substituting  (A-24),  (A-25),  (A-26)  and  (A-27)  into 
(A-22)  and  simplifying,  the  following  expression  is 
obtained: 
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yRG 


=    C 


RG 


4=1  Y    L    fc 
V  K   /      R     Nl 


.8 


8 


(uCp)'2     1 
d  4.8     P2 

i 


T   +r2    -£   ' 
c  q3 


,4-8 


where  it  was  assumed  (uC  ) 

^  p  c 


(viC  ),  at  the  mean 
p  h 


(A-28) 


temperatures  of  the  cold  and  hot  side,  and  all  con- 
stants have  been  written  as  C^r .   The  pressure  term 
in  the  denominator  of  (A-28)  is  the  top  pressure  level 
in  the  cycle. 

Combining  (A-$),  (A-5)  and  (A-16): 


wC 


'x 


P  1-Tl 


X 


.028  (wC  )*8  k  '2 
p  c 

N,8  d  1.8(ptj.4 

c 


d.vSl^i1 


(A-29) 


Divide  equation  (A-28)  by  (A-29)  and  solve  the 
resulting  expression  for  Nd2  : 


wC 


*±  ■  CND  H.5  ^ 


.5 


v  K 


vL.5 


(Pr) 


•3  I 
P 


i+q 


.5 


c      q* 


di 


4.8 


o5 


(A-30) 


Equation  (A-30)  represents  the  cross  sectional  area 
variation  of  the  regenerator  in  terms  of  gas  properties, 
tube  size  and  arrangement,  regenerator  effectiveness, 
flow  rate,  top  pressure  level,  cycle  pressure  ratio,  and 
mean  temperatures  on  the  hot  and  cold  side.   C^,  represents 
the  constant  terms. 
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The  heat  transfer  area  variation  (Nd.L)  is  obtained 
as  a  function  of  the  same  properties  by  solving  (A-29)  for 
LAL,  making  the  appropriate  substitution  for  (Nd2.  )»2 

and  multiplying  this  result  by  (A-30) : 


Nd.L  =  CW,T  wC 
1     NdL 


x 

[  1-T1. 


1.4,   x1.2 

V  K  , 


(Pr) 


64  (\ 
K 


.2 


.8 


.& 


.8 


1.4 


T  +r' 
c 


4.8 


0.4 


(A-31) 


It  is  noted  that  equations  (A-JO)  and  (A-Jl)  simplify 
considerably  for  a  given  gas,  tube  size  and  tube  arrange- 
ment.  The  resulting  expressions  then  become  functions  of 
flow  rate,  regenerator  effectiveness,  top  pressure  level, 
pressure  ratio  and  mean  temperatures  on  the  hot  and  cold 
side.   In  the  last  two  equations,  the  gas  properties  may 
be  evaluated  at  the  mean  temperature  on  the  hot  or  cold 
side  since  Pr  and  the  ratio  of  specific  heats  can  be 
considered  constant  over  a  large  range  of  temperatures 
for  most  gases.   Also  the  change  in  thermal  conductivity 
to  the  0.2  power  would  vary  only  slightly  since  k  under- 
goes little  change  with  limited  temperature  variation. 

III.  Coolers 

Development  of  the  corresponding  cross  sectional 

and  heat  transfer  area  equations  was  done  using  the 

same  analogy  as  for  the  regenerator.   For  the  coolers 
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the  heat  transfer  coefficient  was  based  on  the  outside 
of  the  tubes  (gas  side).   With  the  approximation  that 
the  tube  wall  and  water  film  thermal  resistances  are 
small  compared  to  the  gas  thermal  resistance,  the  over- 
all heat  transfer  coefficient  was  based  solely  on  the 
gas  heat  transfer  coefficient. 

For  the  precooler,  the  cross  sectional  area  vari- 
ation is  then  given  by: 

The  heat  transfer  area  variation  is: 

Similarly  for  the  intercooler,  the  following  results: 

where  R,  K,  k,  and  Pr  are  the  familiar  gas  properties.  In 
the  above  expressions,  the  gas  properties  are  evaluated  at 
the  mean  gas  temperature  in  the  respective  cooler. 

T   =    ^    1   and  T   =   2  P   ?   =   2  P   X 
pc       2  ic       2  2 
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P  and  r  are  again  the  top  pressure  level  and  pressure 
ratio  respectively  of  the  cycle.   For  the  precooler: 

A  Tg  =  T9  -  Tx 

and  A  T  is  log  mean  temperature  difference  across 
the  precooler. 

For  the  intercooler; 


A  Tg  -  T2  "  Tj   -  T2  "  Tl 


andAT  is  the  LMTD  across  the  intercooler, 
m 
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APPENDIX  B 
TURBO -MACHINERY  EQUATIONS  DEVELOPMENT 

I.  General 

The  assumptions  involved  in  the  development  of  the 
rotating  machinery  equations  include: 

1.  Uniform  axial  velocity  over  the  length  of  the 
component. 

2.  All  machines  of  free  vortex  design. 

3.  Equal  work  done  per  turbine  or  compressor  stage 

4.  Static  pressure  and  temperature  equal  to  stag- 
nation pressure  and  temperature. 

5»  Steady  flow. 

Development  of  the  rotating  machinery  equations 
are  in  part  based  on  reference  [6]. 

II.  Compressors 

The  root  stress  in  a  blade  caused  by  centrifugal 
force  can  be  expressed  in  accordance  with  references 
[6]  and  [8]  as: 

(it)  CB-D 


/?  JV        da  (£_ 


g    2 
&o 

where  f  is  the  blade  taper  factor  and  /the  density  of 

blade  material.   Defining  a  stress  parameter  (P)  as 

equal  to  (t~c>)   *  (B-l)  is  written  as: 
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Since  the  annular  flow  area  is: 

A  =   Tidi 
(B-2)  can  be  written  as: 

J\2    A 


(B-3) 


?  = 


4  IT 


2  7ig/ 


2  g   v  d 


rf; 


(B-4) 


where  the  relation  -J\    =  ^y2  was  used  to  obtain  the 
term  on  the  right  in  (B-4-).   From  continuity  and  the 
perfect  gas  relationship: 

AC   =  ^ 
AUa     P 

Using  (B-4),  (B-5)  and  the  relation  RPM  =|°    , 
machinery  rotational  speed  can  be  expressed: 


(B-5) 


RPM  -'  30 


tc  wRT 


2  s0J  Ps" 


17737 


•  5 


(B-6) 


From  (B-6)  it  is  seen  that  RPM  may  be  limited  because  of 
stress  considerations.   For  a  fluid  with  low  sonic  velocity, 
RPM  will  be  limited  by  mach  number  considerations .   For 
RPM  designed  with  mach  number  considerations,  equation 
(B-6)  becomes: 


RPM 


60 


Kg, 


71 


u 


KgQRT 


(i~)    <pca) 


.5 


(B-7) 


after  substituting  (B-4-)  for  P.   All  remaining  turbo 
machinery  relationships  will  be  based  on  stress  considera- 
tions.  For  helium  as  a  working  fluid,  compressor  and 
turbine  design  would  be  influenced  by  stress  considerations 
rather  than  mach  number  limitations. 
(B-3)  can  be  written  as: 

..  I 
d 
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A  =   Ted2 


(B-8) 


(      A   V5 


•=rfj£)  .  (B-9) 

Using  (B-4)  and  (B-5)?  the  pitch  diameter  is  represented 
by: 

d  = 


2  wRT 


2«0p  4 


25 


(B-10) 


For  the  LPC   inlet,    (B-10)   becomes: 

5  /  \     c 


dLPC   =   K 


w"T1r 


a/U' 


i^    -25 
d 


(B-ll) 


where  P  and  r  are  the  top  pressure  level  and  the  pressure 
ratio  respectively  of  the  cycle,  and  K  represents  the 
constant  terms. 

Substituting  (B-ll)  into  (B-8),  the  inlet  annulus 


area  for  the  LPC  is : 


ALPC  "  K 


tfT,r 


O) 


-.5 


(i) 


.5 


1 


(B-12) 


For  most  compressors,  P,   -j  ,  and  C  /  can  be  assigned 
fixed  values  over  a  range  of  operating  conditions.   Then 
pitch  diameter  and  annulus  area  are  dependent  on  the 
volumetric  rate  of  flow  and  cycle  pressure  ratio. 

Using  the  same  procedure,  the  corresponding  equa- 
tions for  the  HPC  are: 

5 


d   -  k/-1 

dHPC  ~  \    P 


(r) 


25 


5  L£\  -25 


wT, 
AHPC  =  K  ~P~   ^r' 


a/U> 


a/U 


n 


O) 


.5  /£ 


(B-13) 


(B-14) 


The  tip  diameter  at  the  inlet  end  of  the  LPC  can  be 
expressed  as 

d 


\  - 


1 


•(I) 


LLPC 


(B-15) 
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with  a  similar  expression  for  the  HPC. 

Applying  Newton's  law  in  the  form  of  angular  momentum, 
along  with  the  assumptions  given  in  section  I  of  this  ap- 
pendix, the  work  done  per  compressor  stage  is: 

T2     C_ 


W 


IT 


s 


■§   (tan  P2  -  tan  ^±)  (B-16) 


where  &o  an<^  ^i  are  ^^-e  angles  as  shown  in  Figure  B-l . 

For  the  LPC,  the  work  is  given  by: 

J  C  T,      c 
WLPC  =   -1^     Ct-5  -  1)  (B-17) 

Dividing  (B-17)  by  (B-16)  the  number  of  stages  is: 
JC  T,(t'5  -  1) 

nlpc=  ri    u-  c <B-18> 

C  SC  i"  ~U  (tan  ^2  "  tan  ^1} 
where  "n    is  the  stage  efficiency.   (B-18)  is  in  the 
form  whereby  mach  number  considerations  might  dictate 
the  number  of  stages  since 

^  -   KRgoTl   -  _JL 
U2      KRg  U"     KRg^M7" 

By  substituting  for  U2  from  equation  (B-4-),  (B-18)  may 
be  put  in  the  form: 

2JC  Tn(T°5-  1)  U  /d), 

NLPC  =  C (B-19) 

T)cnscP  _£  (tan  ^  -  tan  ^) 

In  (B-19),  stress  considerations  now  influence  the  number 
of  stages. 

For  the  HPC  an  expression  similar  to  (B-19)  is 
applicable  since  the  inlet  temperature  to  both  compressors 
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is  equal,  and  only  the  t  ratio  and  the  angles  pp  and- 
P1  may  vary  slightly  from  that  of  the  LPC, 


Compressor 


Turbine 
Figure  B-l 
Stage  Velocity  Diagram 


-77- 


III.  Turbines 

Before  developing  the  corresponding  equations  for 
the  turbines,  it  is  necessary  to  determine  the  pressure 
ratio  across  each  turbine  as  a  function  of  the  overall 
cycle  pressure  ratio.   Since  the  HpT  drives  only  the 
two  compressors,  the  work  of  the  HPT  must  equal  the 
total  compressor  work: 

WHPT  ■  Vp(T6  -  V  -  VpT6  (l  -T^)  <B-20> 

From  equation  (2-2): 

«c  "   2  V1   (t-5.  i)  <B-21) 

Equating  (B-20)  and  (B-21)  and  solving  for  ^Hpm» 

"HPT=        1    -    2T1  ~5     ~    '    1    -    CfT-g    -    I)        (B_22) 

T6"t\ 

where  C  is  a  constant  since  T,  ,  T^,  r\,    and  ti   will  be 
fixed  for  a  given  cycle. 

From  equation  (2-4) : 

(TH.pr)(TLpT)  =  tt  =  Yt  (B-23) 

Substituting  (B-22)  into  (B-23): 

TLfT  =  Yt[1  -  C(t°5  -  1)]  (B-24-) 

At  the  HPT  exit  (LPT  inlet)  the  temperature  can  be 
expressed  as  a  function  of  T^  and  pressure  ratio.  Using 
the  definition  of  turbine  work  and  (B-22): 

T?  =  T6  [1  -  TitC(T'5  ..  i)]  (B-25) 

Following  the  procedure  used  for  the  compressors 
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and  the  preceding  work  in  this  section,  the  pitch  dia- 
meter and  annulus  area  at  the  HPT  exit  can  be  written: 

^  =  " (p[i-c(x-5-i)]»  /  M  W       (B'26) 

2T,  N  x 

where  C  = 


Vt*. 


wT6[i-ntc(T-5-p]      1    (^m5fiy? 


AHPT  =  K    5 {I C     (P)     d7     (B-^) 

HPT     P[1-C(t'5-1)]u      Ca/U       ld/ 

For  the  low  pressure  turbine: 

^-.fr"-'/1-fe"rJ-fe;)-,fff-g    <«»> 


liPT  P  C  /TJ  VK/    Id 

\  a/U/ 

The  tip  diameter  at  the  exhaust  end  of  the  LPT  is 

given  by: 


Dt  - 11  +id/8;  4m 


For  the  HPT,  the  number  of  stages  is 

.  JCpT6  f1-4r")^st 


N^  =      °  v    lHPT  /   ;  bU  (B-31) 


"HPT 


r2      C. 


2-   ^(tan  p2  -  tan  Px) 


w 


ith  (B-22)  substituted  for  Tjrpms  (B-$l)  becomes  : 

H    -   2JCpT6C(T^-l)V8t    /£) 
ft  _a  (tan  ^  "  tan  P])     7 

where  (B-4)  was  used  to  substitute  for  U2  in  (B-31),  C 
is  again  the  term  defined  with  (B=22)o 

For  the  LPT,  the  number  of  stages  become: 

-79- 


NjpT  =    P?i     TLFT  /  t  st  (B-53) 

|-  -§  (tan  P2  -  tan  ^) 


Using  (B-4-),  (B-24),  and  (B-25) 


5  i  m  rv-rCi  rr-r"  5 


N 


P  -§(tan  ^2  -  tan  P1)Yt[1-C(t° -"-I)] 

The  turbo-machinery  equations  developed  in  this  and 

the  preceding  section  appear  to  be  quite  cumbersome.   How= 

ever,  it  should  be  noted  that  many  of  the  factors  (effi- 
'a 


Jt  c 


ciencies^T,  n- ,  P*  Pp  an(^  ^\)   will  be  fixed  by  the  tur- 
bine and  compressor  design,  and  will  be  applicable  to  a 
wide  range  of  operating  conditions,   '■'■'hen  the  only  vari- 
ables left  to  the  marine  engineer's  discretion  are  Tfi, 
T, ,  t  or  r,  Y  and  volumetric  flow  rate. 
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APPENDIX  C 
EVALUATION  OF  CONSTANTS 

I.  General 

The  MGCR  reference  plant  has  a  reactor  rating  of 
49.7  MW  (thermal)  with  a  flow  rate  of  74.6  lb/sec. 
The  coolant  top  and  low  temperatures  in  the  cycle  are 
respectively; 

T6  =  1300°F  and  T±   =  100°F. 
Reactor  inlet  temperature  is  T^  =  791°F.   The  ratio 
of  turbine  pressure  ratio  to  compressor  pressure  ratio 
is  0.86,  giving  a  pressure  drop  of  140  psi  through  the 
cycle.   Compressor  pressure  ratio  is  2 .6  with  the  top 
pressure  of  the  fluid  at  1000  psia.   Turbine  and  com- 
pressor efficiencies  are  respectively; 

tk  =  0.90  and  n   =  .860 
t»  c 

With  the  above  conditions,  the  LPT  produces  20,600  HP., 
and  the  HPT  74-6  auxiliary  horsepower  in  addition  to 
driving  the  compressor  unit.   In  this  work,  the  HPT 
provides  only  sufficient  power  to  overcome  the  com- 
pressor work,  leaving  the  LPT  with  an  output  of  21,500 
HP  for  the  above  conditions. 

The  remaining  temperatures  of  the  reference  cycle 
are  given  below.   The  state  points  correspond  to  those 
indicated  in  Figure  II-l. 
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T1    =  T   =   100°F  Tg  -  852°F 

T2  =  T4  =   239°F  Tq  =  301°F 

T^  =       791°F  Twl  =   70°F 

T6  =      1300°F  Tw2  =   90°F 
T?  =      1008°F 


II.  Regenerator 

The  regenerator  is  composed  of  three  equivalent 
units  arranged  in  series.   Since  no  heavy  internal 
supporting  structures  are  necessary,  the  regenerator 
can  be  constructed  in  a  U-shaped  fashion  to  conserve 
space  when  a  high  regenerator  effectiveness  is  needed. 
This  shape  unit  would  also  simplify  the  thermal  expan- 
sion problems.   The  constants  evaluated  for  the  regenera- 
tor dimensions  and  weight  are  based  on  the  overall  size 
and  weight  of  the  regenerator.   Unit  length  and  weight 
can  then  be  found  by  dividing  by  3.   Each  unit  uses 
3/8  inch  outside  diameter  tubes,  20BWG,  with  a  5/8 
inch  center-to-center  triangular  tube  spacing.   Heat 
transfer  and  cross  section  area  and  effective  tube  length 
are : 

NdL  =  23,340  ft2 

Nd2  =  7.1  ft2 

L   =  81  ft 
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Total  weight  of  the  reference  regenerator  having  an 
effectiveness  of  0.90  is: 
WT  =  37.1  tons 
For  the  above  tube  data,  q  as  given  by  equation  (A-9A) 
in  Appendix  A  is  equal  to  2.06.   With  the  above  infor- 
mation, equation  (2-9)  is  written: 

Nd2  =  KND  w^  j-   (Tc  +  .042l(r)2  T^'5      (C-l) 

with  T  and  T,  defined  by  (2-11)  and  (2-12)  respectively. 
Using  (C-l)  and  the  above  information  KND  is  evaluated 
as: 

%D  =  -91  (0-2) 

Upon  using  (2-10),  KjxDt  becomes: 

KNDL  "  225-5  (C-3) 

Number  of  tubes  is  given  by: 

«■    cross  sectional  area  .     Nd2 /«  ^\ 

£  d'Q  (4  +  1)        ~    l   d!Q  (q  ♦  1) 

where  the  denominator  equals  2.35  x  10  y   for  the  given 
tube  data. 

Effective  tube  length  can  be  found  by  using : 

T    heat  transfer  area      NdL (n   c\ 

L  "    71  N  d  "    71  N  d  ^=>; 

o  o 

Total  regenerator  weight  as  given  by  (2-16)  is: 
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WRQ  =  K-^NdL)  +  K2 

where  NdL  is  equation  (2-10).   Assuming  headers  and  sup- 
port equipment  weigh  3  tons  per  unit,  K~  is  taken  as; 

K2  =  9o0  tons  (C-6) 

Using  (2-16),  (C-3)  soa.6.   the  given  operating  con- 
ditions, K,  is  then  evaluated  as; 

Kx  =  1.21  x  10~5  tons/ft2  (C-7) 

Cross  sectional  and  heat  transfer  area  and  the 
regenerator  weight  can  then  be  written  as; 

.5 


Nd2 


=  .9Kw)(^f  j       f  (Tc  +  .042l(r)2  Th)'5   (C-8) 


,t\        \1.4  ,  .4 

NdL  =  223.5  wM~)      L  Q   (Tc  +  .042l(r)2Th)    (C-9) 

'x      P 


WRG  -  1,21  x  10"5(NdL)  +  9  (C-10) 

with  NdL  given  by  (C-9)» 

III.  Coolers 

Both  coolers  are  similar  in  arrangement  and  con- 
struction.  Each  uses  5/8  inch  outside  diameter  tubes, 
18  BWG,  with  a  7/8  inch  center-to-center  triangular  tube 
spacing.   Weight,  effective  tube  length  and  cross  sec- 
tional and  heat  transfer  area  for  the  precooler  used  in 


-84- 


the  reference  plant  are  respectively; 

W   =  15.6  tons 
pc 

L   =  19.2  ft 

Nd2  =11.6  ft2 

NdL  =  7600  ft2 

Using  equation  (2-17),  the  above  information  and 
the  conditions  stated  in  part  I  of  this  appendix,  KNI) 
is  evaluated  as 

KND  =  1.48  (0-11) 

With  (2-18),  Kjr-jyr  becomes; 

KNI)L  =  290  (0-12) 

Number  of  tubes  and  length  are  given  by  expressions  simi- 
lar to  (C-4)  and  (C-5)  respectively, 

Precooler  weight  is  found  by  using  (2-35)  '• 

.rk6-   /T  V4  (*■]'      4-  Ko         (2-35) 


with  z\T  and  T   given  by  (2-19)  and  (2-20)  respectively, 
g      pc 

AT   is  the  log  mean  temperature  difference. 

Assuming  a  header  weight  of  5  tons; 

K2  =  5  tons  (C-13) 

Using  (2-35),  (C-12),  (C-13)  and  the  given  conditions, 


K,  is  evaluated  as  % 

K1  =  1.13  x  10~5  tons/ft2  (C-14) 

This  value  of  K,  is  of  the  same  order  of  magnitude  as 
that  for  a  steam  condenser  [91  hut  somewhat  smaller. 

Weight,  effective  tube  length  and  cross  sectional 
and  heat  transfer  area  for  the  intercooler  of  the  ref- 
erence design  are; 

W.    =  1$.4  tons 

ic      y 

L  =      14.0  ft 

Nd2      =      13.4  ft2 
NdL     =      6690  ft2 
With   (2-24)   and  the  given  conditions: 

%D  «   3.21  (C-15) 

Given   conditions   applied  to    (2=25)    evaluates  K^rm-    as» 

%DL^77  (0-16) 


Equation  (2-36)  gives  intercooler  weight  ass 

At  \Jo4    .,    .4 

Wic  ■  Klw  hf)        Tic     ^8      +  K2  <2=36) 

\   m/  p 

Kp  for  the  intercooler  is  considered  larger  than  that 
for  the  precooler  because  of  the  higher  pressure  associ- 
ated with  the  intercooler. 

K2  =  6.5  tons  (C-17) 
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Combining  (2-36),  (C-16),  (C-17)  and  given  conditions 

K±   =  1.03  x  10~5  tons/ft2  (C-18) 

Summarizing,  the  equations  for  the  precooler  and 
intercooler  ares 

Precooler; 

/AT  y5 

Nd2  =  1.4-8  wf  -=£)   t   °5  £  (C-19) 


NdL  =  290  w^"   V   (pJ         (C^0) 

W   =  I.13  x  10~5  (NdL)  +  5  (C-21) 

pc 


with  NdL  given  by  (C-2  0). 
Intercooler; 

Nd2  =  3.21  wH^)   T.c'5  ^        (C~22) 

NdL   -   477  wLL.1  T  *_  (C=23) 

UT    /  1C         p. 8 

\*    m  * 

W.      =    L03  x   1(T5    (NdL)    +   6.5  (C-24) 

J- '- 

with  NdL  given  by  (C-23). 

IV o  Reactor 

Core  radius  cubed  for  the  reference  reactor  was 
developed  in  sub~section  C  of  Section  II  and  can  be 
written  from  equation  (2=52)  ass 


-a 


n 


R3  =  KR  2-g   (T6  -  T^)1'4  (C-25) 

With  the  core  diameter  of  the  MGCR  reactor  equal  to  6.87 
ft,  and  for  the  given  operating  conditions,  KR  is  evaluated 
from  (C-25)  as: 

KR  =  .0176  (C-26) 

The  weight  of  the  reference  ractor  is  approximately 
173  tons.   Combining  (2-51)  and  (2-53): 

WR  =  (\   5.81  R3  (C-27) 

Substituting  (C-25)  and  (C-26)  into  (C-27)  and 
applying  the  given  conditions,  reactor  density  is  evalu- 
ated ass 

(°R   -  .92  tons/(ft3  of  reactor  core)  (C-28) 

Considering  the  reactor  vessel  (exclusive  of 
shielding),  f°R  =  .105  tons/ft3.   Substituting  (C-25), 

(C-26)  and  (C-28)  into  (C-27) ,  reactor  weight  becomes  § 


WR  =  .094-5  (T6  -  T5)le4  -2g  (0-29) 


Core  radius  cubed  is  given  by; 

R3  -  .0176  ^q   (T6^T5)le4  (C-30) 

Core  height  can  then  be  found  using  (2-51) • 

H  =  1.85R  (2-51) 
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Reactor  vessel  radius  (Ry)  and  height  Hy  are  given 
by  (2-56)  and  (2-57) : 

Ry  =  1.65  (R)  (2-56) 

Hv  =  5.51  (H)  (2-57) 

V.  Turbo  Machinery 

In  developing  the  turbo  machinery  equations,  a 
large  number  of  variables  are  present.   Fortunately  most 
of  these  variables  are  a  function  of  the  machine  design 
and  not  a  strong  function  of  the  cycle  operating  condi- 
tions.  Therefore  most  similar  machine  variables  will 
have  values  that  are  near  those  indicated  below  (for  the 
reference  plant)  for  a  considerable  range  of  operating 
conditions. 

Both  machines : 

P  ^>    10,000  -  12,000  ft 
Ca/U  ^  1/2 

Compressor      Turbine 
t  (low  pressure  end)  £&  .11  .20 

efficiency  C&  .86  .90 

Weight,  machine  length  and  tip  diameter  of  the 
rotating  machinery  are  considered  to  be  of  significance 
in  determining  the  overall  plant  size  and  weight.   There- 
fore the  constants  associated  with  each  of  these  three 
developed  equations  (weight,  length  and  tip  diameter)  were 
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evaluated  for  each  rotating  machine. 

The  following  data  apply  to  the  machinery  used  in 
the  reference  design: 


LPC 

HPC 

HPT 

LPT 

L   (ft) 

7.95 

10.25 

3.98 

2.84 

W   (tons) 

4.8 

3.92 

3.42 

3.33 

D   (inches) 

18.1 

14.3 

20.2 

25.0 

The  applicable  length  (L) ,  weight  (W)  and  tip  dia- 
meter (D)  equations  from  sub-section  D  of  Section  II  were 
used  in  conjunction  with  the  above  given  geometry  data 
and  the  operating  conditions  to  evaluate  the  constants. 
The  following  values  were  obtained  as  constants: 

LPC      HPC       HPT       LPT 
Ko     (ft)         37.0     47.6      18.5      27.2 

K  Lbf  sec       113.0    151.0     156.0     216.0 


W      .   2 

in 


5 


Kd  f  GLbmbf  > 

For  the  LPC,  the  working  equations  are; 

L  =  37.0(t°5  -  1)  (C-51) 

W  =  115.0  |F  (t°5  -  1)  (C-32) 

D  -  40.7  (f)'3  (C-33) 

For  the  HPC: 

L  =  47.6(t°5  -  1)  (C-34) 
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W  =    151.0  §  r°5    (t*5  -   1) 


P 


/|)°5     r' 


For   the  HPT; 

L  =   18.5    (t'5  -   1) 


W  .  156  1  (^5-i)[i-y(^3-i)3 

[1   .   c(tc5   -   l)]u 


2V5   /[I-tj^CCt"5-!)]   V 


(c- 

=35) 

(C- 

-36) 

(c- 

=37) 

(c 

=38) 

D  .  63.4(f)'5  fCi-ntcCT^-i;j     -  (c„39) 


[1   -   C(t°-1)]u 
2T 


1  K 

where   C   =   a  _   _      and  u  = 


T6Vc  K=1 


For   the   LPT 
27o2    [l-TitC(To5-l)][YT(l-C[T'5-l])~l] 
Yt    [1-C(t°5-1)] 


l  = £— — ■= — : : (c-40) 


W  =   216.0      EE 


[l-Tit(l^)][l~TitC(To5»l)][YT(l«C[To5-]])-l]      (C-/+1) 


Yt[1   -   C(t°  -  1)] 


/jblV      n   -  «     n    .     i-^i-5 


D  =   65.3    (-f-l        [1   -  r,t    (1   -     £-.)]•?  (C~42) 
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APPENDIX  D 
SAMPLE  CALCULATIONS 

I.  General 

In  order  to  indicate  the  effect  of  certain  para- 
meters as  they  individually  influence  plant  efficiency, 
weight  and  space,  it  is  desirable  to  fix  as  many  other 
variables  as  possible.   As  stated  in  Section  I,  the 
top  and  low  coolant  temperatures,  reactor  inlet  tempera- 
ture, turbo-machinery  efficiencies  and  pressure  losses 
for  the  cycle  will  generally  be  fixed  for  a  power  plant 
of  this  type.   From  consideration  of  the  developed  equa- 
tions, the  remaining  parameters  effecting  the  cycle  are 
pressure  (P) ,  flow  rate  (w) ,  pressure  ratio  (r) ,  and 
regenerator  effectiveness  (r\    ).   However,  as  noted  in 
Section  I,  regenerator  effectiveness  is  specified  once 
a  pressure  ratio  is  selected  for  this  type  cycle,, 

Solution  of  all  equations  was  performed  for  a 
LPT  (net  power)  output  of  21,500  HP  with  the  below  fixed 
conditions s 

T.   =  1760°R  r\.       =  .90 


6  -  *'*~  -  '"t 

\      -      560°R  i\c 


h 

=  1251°R 

p 

=   1000  psia 

t  1 

70°F 

r 

;w2 


=  .86  or  Y  =  .94-1 


t„o   =  90°F 
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Another  set  of  calculations  were  made  with  the 
reactor  inlet  temperature  lowered  to  1160  R,  with  all 
other  parameters  above  held  constant. 

Equations  were  evaluated  for  pressure  ratios  rang- 
ing from  1,3  to  a  ratio  that  placed  the  LPT  exhaust 
temperature  (Tg)  at  the  same  level  as  the  reactor  in- 
let temperature  (T(-).   This  condition  would  require  a 
regenerator  effectiveness  ("n  )  of  100  yo.  AnY  pressure 
ratio  larger  than  this  value  would  present  an  impossible 
cycle  since  r\     would  have  to  be  greater  than  100  /o. 

For  equations  where  w  and  P  appear,  w  is  in 
lb/sec  and  P  is  in  p  sia.   All  temperatures  are  in 
degrees  absolute  (  R) . 

II .  Sample  Calculations 

Sample  calculations  are  indicated  for  all  equations 
using  a  pressure  ratio  of  2.4  and  the  fixed  conditions 
given  in  Section  I  of  this  appendix. 

A.  Flow  rate 

Flow  rate  is  found  using  equation  (2-5) ' 

T  C 
net  power  -  w  [t»t  +  T&C  (1-  ^)-  2  ,  X  p  (t°5  -  l)] 

c 

^'^^^.90(i76o)i.a^(^ntngJ-aft6o^^Vi.i9-i)3 

w  =  81.1  lb/sec 
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w    = 


%ili55002    .    i3.6ib/fHT>.m, 


21,500 


B.    Efficiency 

Using    (2-8),    efficiency  becomes 

i  =  °90  ^t§  C1  -  nfeiaaj  -  t^1-1*  - 1} 

1760  -  1251 
560 

i\   =   0.294 

Co  Cycle  Temperatures 

T4  is  found  using  the  following  equation  obtained 
from  compressor  works 

T4  =  Ti  [1  +    vT  (t°5  "  1)] 


T4  =  T2  =  560  [1  +  -|g  (1.19  -  1)]  =  686°R 


Tg  is  given  by: 


T8  "  T6  [1  -  V1 


Yt 


)] 


=  1760 


1  -  .90  1  - 


(.941)1.42J 


T8  =   1562UR 


Tq  is  found  from  a  heat  balance  around  the  regnerator; 


■94- 


T5  -  T4  =  T8  -  T9 


Tq  =  1362  -  1251  +  686 


t9  =  797°R 


D,  Regenerator  effectiveness 

ti   is  given  by  (A-4-)  in  Appendix  A; 


T5  -  \ 


t,      _      1251   -  686 
^x  "      1362  -  686 


tjx  =      0.856 


E.  Regenerator  weight  and  dimensions 

T  =  T5  +  T4  =   1251  ±   686  =   968°R 
c       2  2 

Th  .  !s_p2  .   1?62  +  797  .   108qOR 
Putting  numbers  in  (C-8),  (C-9)  anci  (C-10)  respectively t 

Nd2  =  .91(81.1)  (j^fgg)  '  10OO  (968  +  •04-21(2.471080)°' 
Nd2  =  5.85  ft2 

NdL  =  223.5  (81.1)  (XI7836)  0(1000)°8  (968+. 04-21(2,47 1080) 
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NdL  =  12,210  ft2 

WRG  =  1.21  x  10~5  (12,210)  +  9 


WRG  =  23»8  tons 


Using  (C-4-)  and  (C-5)»  the  number  of  tubes  and 
length  are  respectively? 


N  =  5^2 m        2^90  tubes 


2.35xlO~5 


L  =   12*210 ^   =   50  ft 

8/  12 


(3.1^)(2490)f|)  i- 


l-    5a- 


3  =   16.67  ft/unit 


F.  Precooler  weight  and  dimensions 

T   .  !2^_!i  .  797  +  ^60  .  ^o^ 

pc       d  Z  — ' — — — 


A!  =  T9  -  ^  =   797  -  560  =  237 


A  T  =   T9  I  Tw2  '  (T1  "  Twl}  =  797  -  550  -  (560  -  550) 

J   n   T9  "  Tw2  ^   2i£Z 

-^  n   T     T  ^  n   30 

1    wl  y 


/\   T   =  105°F 

m   —— ■* — 


>96= 


Using  equations  (C-19),  (C~20)  and  (C-21) 
Nd2  '  *5  ™'5 


=  1.48  (81.1)  (f^)  "   (678) 


1000 
N<1'  =  11.58  ft' 

NdL  .  290  (81.1)(^2)ls4  (678)-4  (fs^o)'8 
NdL  =  8220  ft2 
Wpc=  1.13  x  10"3  (8220)   +  5 


wpc  =  14 > 5  tons 


N   =  11.^8 
.00461 


N   =  2470  tubes 
8220 


L 


x   (2470)  (|)  ±2 
L   =  20.4  ft 

G.  Intercooler  weight  and  dimensions 


T   +  T 

m   =    2   1    686  +  560     c-ozOD 

xic      2     =  ? —   =  bed*   R 


A  T  =   T2  -  Tx  =  686  -  560   =  126 


-97- 


A  T  .  686  -  ^0  -  (^60  -  ^0) 

^  m  "       n  1^6 

A  Tm  =  70°F 

Using  (C-22)  through  (C-24-); 

Nd2  =  3.21(81.1)  (if§)°5  (623)e5  ^fg^ 

Nd2  =  15.5  ft2 


NdL  =  ±77   (81.1)  (±§5)  *   (623)  *4  (2,4?'o 

w°'  (1000)  °8 


NdL  =  6590  ft2 


W  .   =  1.03  x  10~3  (6590)  +  6.5 


1C 


W  . 


ic  =  13Q  tons 


N 


=  ^g|61  =   22iO_tubes 


L   =  <^Q- 


71(294-0)  |  ^2 

L   =  15.7  ft 

H.  Reactor  weight  and  dimensions 
Using  (C-29)  and  (C=30): 


Wp    o094-5  (1760  -  1251)1'4   81.1 

~98- 


R  8 

(1000) °° 


WR  =  188  tons 


R3  =  o0176  (81al)   8   (1760  -  1251) 


1.4 

(1000) 


R3  =  35^04 

R  =  5.28  ft  (core  radius) 
From  (2-51)  core  height  is; 

H  =  1.85  (5.28) 

H  =  6.06  ft   (core  height) 

Reactor  vessel  radius  and  height  are  obtained  from 
(2-56)  and  (2-57)  respectively; 

Rv  =  1.65  (3.28) 

Ry  =  -2»  ^?   ft 

Hv  =  3.31  (6.06) 
Hy  =  20  ft 

I.  Rotating  machinery 

Using  equations  (C-31)  through  (C-42),  the  lengthy 
weight  and  tip  diameter  are  given  respectively  for  each 
unit ; 

LPC? 

L  =  37  (1.19-D 

L  =  7.0  ft 

w  =  113  (8i6oo(2^)  (i-19-i) 
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W  =   4.2   tons 

D  .  4.0?     (sL>U2Ja\  °5 
V  1000 

D  =  18  inches 

HPC: 

L  =  47.6  (1.19-1) 
L  =  9.1  ft. 

W  =  151   f5§J  (2.4)'5  (1.19-D 
W  =  5.61  tons 

D  ■  tt-x  (low")'5  ^2-*>*25 

D  =  14.5  inches 

HPT: 

L  =  18.5  (1.19  -  1) 
L  =  5.52  ft. 


w  .  156  81^1  (1.19-l)[l-»9(-822)(l. 19-1)1 
1000    [1-. 822(1. 19-1)]2'5 


W  =  5.2  tons 

D   -   6z  4  fffiLxlY5   /[l-.9(-822)(1.19-l)]N 
'   fe5.4  ^1Q00y        ^[1.o822(l#19.l)]2.5 

D   =   20.7   inches 

-100-= 


LPT 


L  =  27.2 


[l-.  9(.  822)  (1. 19-1)  ][.  94-1(1. 42)(1-.822[1. 19-1])] 


941(1.42) [1-. 822(1. 19-1)] 


L  =  2.8  ft 


W  _  pi  c.   81 . 1^2 .4^ 

-  ^J-o    1000 


F"9(  1"7941(T742)  }  ]  [l--9(.  822)  (1. 19-1)  ]C 941(1. 42)  (1-.822DLJ.9-Ij 


. 941(1. 42) (1-. 822 [1.19-1]) 


W  =  5.2  tons 


D  =  65o3 


(81.1)2.4 
1000 


.5 


1  - 


.941(1.42) 


D  =  25.4  inches 
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